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ABSTRACT 

Carmichael  C1H  has  proposed  a  direct  injection, 
temperature  ignition  <TI),  reciprocating  engine  design  which 
divides  the  conventional  -four  stroke  cycle  functions  between 
a  compression  cylinder  and  an  expansion  cylinder, 
interconnected  by  a  static  regenerative  heat  exchanger.  The 
prediction  of  heat  transfer  loss  in  the  expansion  cylinder  is 
required  to  predict  the  performance  of  the  new  cycle. 

Woschni 's  correlations  C9U ,  developed  for  direct 
injection,  quiescent,  four  stroke,  compression  ignition  (CI) 
engines,  are  used  to  predict  uncorrected  expansion  cylinder 
heat  transfer  loss.  Laws  of  similarity  allow,  in  principle, 
extrapolation  outside  the  experimental  data  range  for  forced 
convection  heat  transfer  only.  General  applicability  of 
Woschni 's  correlations  is  established  by  examining  the 
expansion  cylinder  quiescent  condition,  using  predicted  mean 
gas  velocity  magnitudes  and  fully  developed,  turbulent  flow, 
forced  convection  models.  Predicted  radiant  peak  flame 
emissivity  and  temperature  are  found  and  compared  with 
nominal  peak  values  observed  in  conventional  CI  engine  cycles 
to  establish  a  minimum  radiation  flux  reduction  ratio.  A 
mean  reduction  ratio  for  forced  convection  flux  is  also 
found.  A  corrected,  conservative  prediction  of  heat  transfer 
loss  in  the  expansion  cylinder  is  then  established. 
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NOMENCLATURE 

<-> 
heat  transfer  surface  area     Zft^l 

»-> 
cylinder  heat  transfer  sur-face  area  Cft"1"] 

regenerator  port  exit  cross-sectional  area 
C-ft2] 

crank  radius  length  Cft3 

bore  diameter   Cft3 

-fuel  consumption  rate     Clbm/hrU 

Boltzmann  number;  de-fined  by  equation  (61) 

■friction  coefficient;  defined  by  equation  (11) 

specific  heat  at  constant  pressure 
CBTU/lbm-oR3 

specific  heat  at  constant  volume 
CBTU/lbm-°R3 

C-,   Woschni  correlation  constants;  defined  by 
equations  18a  —  lSe 

c  speed  of  light  =  2.997902  x  1010  cm/sec 

cb         Boltzmann 's  constant  =  1.38026  x"  10-16  erg/°K 

c  Planck's  constant  =  6.623777  x  10     erg-sec 

cs         soot  formation  rate  coefficient 
C  lbm/sec-f  f~'-atmD 

D  diameter   [ft: 

E         "radiant  emissive  power  CBTU/hr-f t^  or 
BTU/sec-ft2: 

E.         radiant  monochromatic  emissive  power 
C ergs/sec-cm  -\ (cm) 3 


Al 

phanumer i  c 
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Be 

Bo 

Cf 

'   CP 

Cv   ' 

Cl  ?< 

C 

Es         soot  combustion  activation  energy 
CBTU/moledbm)  1 


F  radiant  geometric  view  factor 

h  heat  transfer  coefficient  CBTU/hr-f t2-°R] 


hQ         heat  trans-fer  cae-f-f  icient  across  pipe^ 

diameter  D;  turbulent  -flow  CBTU/hr  — F  t2-0  R3 

h  total  cylinder  heat  transfer  coef-f  icient; 

turbulent  -flow  CBTU/hr-f t2-°R] 

h-.  heat  trans-fer  coefficient  across  length  1  of 

flat  plate;  turbulent  flow  CBTU/hr-f t2-°R 3 

h  local  heat  transfer  coefficient  at  distance  x 

from  leading  edge  of  flat  plate;  turbulent 
flow  CBTU/hr-f t^-oR] 

j  Colburn  factor;  defined  by  equation  (9) 

K  thermal  conductivity  CBTU/hr-f t-° R3 

kr         apparent   flame  absorption   coefficient   per 
unit  length  Cft   3 

k  apparent  gas  absorption  coefficient  per  unit 

9         length  Cft-13 

k  apparent   flame  absorption   coefficient   per 

unit  length-unit  pressure  Cft-atm3 

k^         monochromatic  absorption  coefficient  per  unit 
length  Cft-AlT1 

L  piston  stoke  Cft 3 

l_£  cylinder  gas  volume  height  Cft3 

Z  flame  luminosity 

1  flat  plate  length  Cft3 

lc  Re  characteristic  length  Cft3 

lp  effective  radiant  mean  beam  length  CftD 

1  connecting  rod  length  CftD 

m  charge  mass  Clbm3 

mc  charge  mass  at  ignition  start  Clbm3 

• 

m  charge  mass  flow  rate  Clbm/sec3 

Nppjvj  cylinder  cycle  frequency  CRPM3 

Nu  Nusselt  number;  defined  by  equation  (3) 

P  cylinder  pressure  CatmH 
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P  cylinder  pressure  at  ignition  start  Catmll 

P  cylinder  motoring  pressure  CatmD 

Pr         Prandtl  number;  de-fined  by  equation  (5) 

P  unburnt   -fuel  partial  pressure  in   the   soot 

formation  zone  CatmJ 

Q  heat  transfer  CBTUII 

Q  heat  transfer  rate  CBTU/sec] 

q  heat  transfer  flux  CBTU/sec-f t^D 

q.  monochromatic  radiant  heat  transfer  flux 

CBTU/sec-f  t2-,\] 

R  crank  ratio  =  1  /a 

R  universal  gas  constant  CBTU/mol e ( 1 bm) — ° R 1 

R  specific  gas  constant  =  53.90  lbf -f t/lbm-°R 

Re  Reynolds  number;  defined  by  equation  (4) 

Ren  Reynolds  number  for  pipe  flow;  1   =  D 

Re         local   Reynolds   number  for   flat   plate   at 
distance  x;  1   =  x 

r_         cylinder  compression  ratio 

5  piston  speed  Cft/secH 

6  soot  concentration  Clbm/ft'"'3 

SFC        specific  fuel  consumption  Clbm  fuel/HP-hr3 

St         Stanton  number;  defined  by  equation  (10) 

s  straight  line  distance  between  cylinder  crank 

shaft  and  piston  pin  CftU 

T  temperature  C°R] 

Tg         cylinder  bulk  gas  temperature  L°R1 

Tgc        cylinder   bulk  gas  temperature   at   ignition 
start  C°R3 

T  cylinder  motoring  bulk  gas  temperature  C°R3 


T  local  temperature  of  the  soot  formation  zone 

C°R3 

Tw         heat  transfer  surface  wall  temperature  [°R] 

t  unit  time  Csec3 

V  cylinder  gas  volume  Cft'~'3 

Vc  cylinder  clearance  volume  at  piston  TDC  Cft'~'3 

V  cylinder  volume  at  ignition  start  (ft"') 

V  volume  of  the  soot  formation  zone  Cft""'3 

U  bulk  fluid  (gas)  velocity  Cft/sec3 

Ui         impulse  charge  velocity  exiting  regenerator 
port  Cft/secD 

Uexp       bulk   gas  velocity  at  the  cylinder  walls   in 
the  expansion  stroke  Cft/sec3 

Uexh       bulk  gas  velocity  at  the  cylinder   walls   in 
the  exhaust  stroke  Cft/sec3 

a>  angular  gas  velocity  at  radius  r  Cft/secD 

X  radiant  path  length  Cft3 

Symbol s 

$  charge  equivalence  ratio 

$  unburnt  fuel  equivalence  ratio  in  the  soot 

formation  zone 

©  crank  angle  C°  degrees] 

&  forced  flow  boundary  layer  Cft3 

X  electromagnetic  wavelength  Ccm  or  ,hjti3 

)*.  dynamic  viscosity  CI  bm/f  t-sec  3 

/>  density  Clbm/f  t'~'3 

•v^.         kinematic  viscosity  (momentum  diffusivity)  = 
>*.//>  C-ft^/sec] 

«  thermal  diffusivity  =  K/C  p  Cft2/sec3 

H"  ratio  of  specific  heats  =  C  /C 
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T  apparent  shear  stress  acting  on  the  -fluid  at 

the  wall  Clb-f/f t2H 

£  emissivity  coef  f  i  ci  ent ;  grey  body 

£  apparent  -flame  emissivity  coefficient;  grey 

body 

£r         flame  emissivity  coefficient;  grey  body 

£  gas  emissivity  coefficient;  grey  body 

£  emissivity  coefficient  of  infinitely  thick 

radiating  flame;  grey  body 

£  wall  emissivity  coefficient;  grey  body 

T  transmi ttance 

i»  absorptance 

!»,  monochromatic  absorptance 

^  reflectance 

^  wall  reflectance 

—9 
a-  Stef an-Bol tzmann  constant  =  1.713  x  10 


o 


BTU/hr-  ft'-°R4 


A  incremental  change  or  difference 

A  heat     specific  heat  release  per  cycle 
CBTU/lbm-charge] 

ft  heat  transfer  loss/heat  release  per  cycle 

V  nominal  Qp-./Qj  per  cycle  in  conventional 

four  stroke,  high  rc ,  direct  injection  CI 
engi  nes 
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Subscripts 

cp  combustion  products 

FC  -Forced  convection  heat  transfer 

NCE  New  Cycle  Engine 
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T  total  heat  transfer 
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6  crank  angle  degrees 
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CHAPTER  ONE 
INTRODUCTION 

1.     SCOPE 

Conventional  compression  ignition  (CI)  reciprocating 
engines  have  -found  significant  application  in  commercial, 
industrial,  and  in  land  and  marine  transportation  systems. 
Recent  emphasis  on  -fuel  costs  and  availability,  in 
conjunction  with  improved  materials  engineering,  has  prompted 
consideration  of  design  alternatives  which  offer  increased 
thermal  efficiency  and  reduction  in  operating  costs.  A  New 
Cycle  Engine  design  has  been  proposed  by  Carmichael  L12  in 
which  the  functions  of  a  conventional  four  stroke 
reciprocating  cylinder  are  divided  between  a  compression 
cylinder  for  charge  induction  and  compression,  an  expansion 
cylinder  for  charge  combustion-expansion  and  exhaust, 
interconnected  by  a  static  regenerative  heat  exchanger. 

The  attractiveness  of  such  an  engine  is  possible 
because  improved  fabrication  techniques  presently  allow 
construction  of  strong,  thermal  resistant,  extremely  thin- 
walled  ceramic  matrices  possessing  high  heat  transfer 
properties.  A  cyclic  interchange  of  compressed  charge  into 
and  expired  exhaust  gas  out  of  the  expansion  cylinder  via  the 
regenerator  establishes  an  effective  medium  for  heat  transfer 
and  allows  use  of  exhaust  heat  otherwise  lost  to  the  system. 
The  temperature  of  charge  air  heated  by  adiabatic  compression 
is   boosted   from  the  regenerator   heat   exchange,   a   direct 
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consequence  of  which  is  significant  reduction  in  cylinder 
compression  ratio,  r  ,  required  to  achieve  optimum 
efficiency.  A  computer-generated  thermodynamic  model  of  this 
engine  in  preliminary  design  indicates  that  improved  thermal 
efficiency  may  be  achieved.  An  important  element  in  the 
cycle  analysis  is  accurate  estimation  of  heat  transferred 
from  the  gas  medium  in  the  expansion  cylinder  to  its 
enclosing  surface  area,  and  ultimately  to  the  coolant.  The 
magnitude  of  this  heat  transfer  reduces  both  heat 
availability  to  the  regenerator  and  the  thermal  efficiency. 

The  purpose  of  this  thesis  is  as  follows: 

1.  To  evaluate  the  applicability  of  existing 
heat  transfer  correlations  developed  to 
predict  total  mean  instantaneous  heat 
transfer  coefficient,  h  ,  in  four  stroke 
direct  injection,  quiescent,  CI  engines 
for  prediction  of  total  heat  transfer  in 
the  New  Cycle  Engine  expansion  cylinder. 

2.  To  evaluate  and  compare  the  magnitude  of 
the  reduction  ratio  in  both  the  forced 
convective  and  the  radiant  heat  transfer 
from  the  high  peak  combustion  pressure 
conditions  inherent  in  high  r  ,  CI 
engines  to  low  peak  combustion  pressure 
conditions  inherent  in  the  low  r  ,  New 
Cycle  Engine  design. 
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3.  To  use  the  reduction  ratios  to  correct 
existing  heat  transfer  correlations,  1 f 
applicable,  and  establish  a  correction 
allowing  better  prediction  o-F  heat 
transfer  loss  in  the  New  Cycle  Engine 
expansion  cylinder. 
2.     NEw1  CYCLE  ENGINE  DESIGN  BACKGROUND 

A  vertical  cross-section  of  the  New  Cycle  Engine  design 
is  provided  in  Figure  1-1.  For  simplicity,  each  cylinder  is 
considered  of  simple  geometry,  with  straight  forward  walls 
and  flat  cylinder  head/piston  crown  configuration;  bore 
diameter,  B,  equals  stroke  length,  L.  Gas  velocities  are 
boosted  by  piston  speed,  induction  suction,  charge  impulse 
exiting  the  regenerator  port, %and  combustion  turbulence  only; 
quiescent  flow  should  predominate  with  no  deliberate  design 
feature  incorporated  to  accelerate  any  rotational  velocity 
incurred  by  off-set  regenerator  port  impulse.  The  cylinder 
crank  shafts  will  be  mechanically  linked  to  enable  the 
compression  piston  to  cycle  at  twice  the  speed  of  the 
expansion  piston,  alternating  between  an  induction- 
compression  and  a  free  motoring  stroke.  Proper  crank  angle 
phase  shift  will  exist  to  sequence  the  compression  piston 
motion  at  the  start  of  expansion  piston  down-stroke  to 
optimize  charge  transfer. 

Correct   value   timing  is  imperative  and   is   sequenced 
through  one  charge  cycle  as  follows: 


1.  Compression  cylinder  induction  stroke: 
induction  inlet  valve  1  opens; 
regenerator  transfer  valve  2  is 
cl osed . 

2.  Compression  cylinder  compression 
stroke:  valve  1  closes;  valve  2 
remains  closed. 

3.  Expansion  cylinder  combustion- 
expansion  stroke:  regenerator  exhaust 
value  3  closes;  valve  2  opens  inducing 
compressed  charge  impulse  through  the 
regenerator   as   the  expansion   piston 

.       starts  the  down-stroke;   valve  2   then 
cl oses. 

4.  Expansion  cylinder  exhaust  stroke: 
valve  3  opens,  regenerating  heat 
exchange  capacity  during  exhaust  gas 
expiration.  Valve  2  remains  closed; 
valve  1  is  cycling  to  support  the  next 
charge  cycle  as  expansion  piston 
reaches  top-dead-center  (TDC) . 

A  vertical  cross-section  schematic  of  the  expansion 
cylinder  and  piston  linkage  showing  principle  dimensions  is 
provided  in  Figure  1-2.  A  derivation  of  geometric 
relationships  required  to  calculate  instantaneous  expansion 
cylinder  heat  transfer  surface  area,  ^HT'  volume,  V,  and 
piston  speed,   Sp ,   is  enclosed  as  Appendix  B.   Instantaneous 
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expansion  cylinder  pressure,  P,  mean  gas  temperature,  Tg ,  and 
volume,  V,  at  5°  crank  angle  increments  -From  the  computer- 
generated  thermodynamic  model  is  inclosed  as  Appendix  A, 
Table  A-l  C13;  associated  P-V  and  Tg-v"  diagrams  -for  the 
entire  cycle  are  illustrated  in  Figures  1-3  and  1-4 
respectively.  Other  principle  design  and  operating 
parameters  -for  the  expansion  cylinder  required  to  support  the 
heat  transfer  calculations  are  listed  in  Appendix  A. 
Appendices  A  and  B  provide  the  data  point  basis  -for  the  heat 
transfer  calculations.  A  consideration  of  piston  ring  land 
AreA,  cylinder  sleeve  clearances,  and  clearance  volume,  Vc , 
at  TDC  is  neglected.  All  the  area  and  volumetric 
calculations  Are  based  on  simple  cylinder  geometry.  Within 
the  accuracy  of  reciprocating  engine  heat  transfer 
correlations  in  general,  errors  associated  with  these 
assumptions  Are    considered  to  be  insignificant. 
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CHAPTER  TNQ 

FORCED  CONVECTION  HEAT  TRANSFER  IN  THE 
NEW  CYCLE  ENGINE  EXPANSION  CYLINDER 

1.     THEORETICAL  BACKGROUND 

When  a  temperature  difference  exists  between  a  solid 
surface  and  a  moving  fluid,  thermal  energy  is  transferred  by 
the  joint  process  of  conduction  and  convection.  Shear 
stresses  acting  on  the  fluid  will  slow  the  fluid  immediately 
adjacent  to  the  surface,  forming  a  thin  boundary  layer  of 
thickness  & .  The  heat  transferred  between  the  bulk  fluid  and 
the  body  is  transported  by  conduction  in  the  boundary  layer 
as  a  function  of  fluid  movement  source  (natural  or  forced) , 
flow  character  (laminar  or  turbulent),  physical  properties  of 
the  fluid,  and  body  geometry  and  roughness.  The  heat 
transfer  rate,  Q,  by  conduction  in  any  direction  x  at  any 
instant  is  determined  using  Fourier's  Law  by  a  heat 
conduction  coefficient,  K,  and  the  temperature  gradient: 

Q  =  -KAVT  (1) 

"Ox 

In  laminar  flow,  fluid  particles  flowing  parallel  to 
the  body  surface  conduct  heat  in  a  normal  direction.  At  a 
certain  transition  velocity,  turbulence  is  created  where 
particles  have  a  velocity  component  normal  to  the  bulk  stream 
direction;  the  intense  mixing  of  individual  layers  promotes 
rapid  transportation  of  heat  to/from  the  surface.  Even  in 
turbulent  flow,   however,  a  thin  boundary  layer  exists  due  to 


■fluid  viscosity  across  which  laminar  -flow  and  heat  transfer 
characteristics  are  maintained.  In  turbulent  systems,  the 
heat  transferred  between  a  body  surface  and  moving  fluid 
remains  limited  by  boundary  layer  thickness,  o*  and  fluid 
thermal  conductivity,  K. 

An  ideal  linear  velocity  and  temperature  gradient 
across  the  boundary  layer  between  a  gas  medium  and  a  body 
surface  allows  approximation  of  the  temperature  gradient  as 
<Tg-Tw)/&;  substitution  into  equation  (1)  yields  the 
Newtonian  expression  for  steady  flow  conditions: 

Q  =  -KAVT 
"dx 

=  -KACTg-Tw] 
6 

=  -K(Tg-Tw) 


where  letting  h  sb  — K  : 

6 


Q  =  hA(Tg-Tw)  (2) 


For  real  boundary  layers  of  varying  velocity  and  temperature 
gradients,  h  may  be  determined  in  some  cases  by  integration 
of  the  differential  equations  describing  the  heat  transfer 
process.  For  more  general  and  practical  applications, 
however,  laws  of  similarity  may  be  used  to  predict  the  value 
of  h.  According  to  the  laws  of  similarity,  it  is  sufficient 
to  determine  the  characteristic  relations  between  independent 
dimensionless  groups  obtained,  recognizing  that  any  such 
relations  will  apply  to  all  similar  heat  transfer   processes. 
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From   the   Buckingham   pi  theorem  which   allows  reduction  of 

initial   variables   to   minimum   dimensi onl ess  groups,    the 
following   di mensionl ess  numbers  ars    used  to  describe   forced 

convection  heat  transfer  in  a  fluid  medium  to  a  solid  body: 

1.  Nusselt  #,  Nu  =  hlc/K  <3> 

2.  Reynolds  #,  Re  =  Ulc/T/k  <4> 

3.  Prandtl  #,  Pr    =  v,./«  =  ^C„/K  (5) 

K  p 

where:      v^   =   >t/p   =   momentum   diffusivity 

(kinematic   viscosity)    (6) 

«  =  K/C  p  =  thermal  diffusivity  <7) 

,n.  =  dynamic  viscosity 

p    =  density 

C   =  specific  heat  at  constant  pressure 

1   =  Re  characteristic  length 
and  in  general  form: 

Ru  =  function  of  Re,  Pr  =  f  (Re,Pr)       (3) 

In  fully  developed  turbulent  flow  systems,  analysis  of  heat 

transfer   processes  is  enhanced  by  a  simple  analogy  between   heat 

and  momentum  transfer  which  allows  prediction  of  h  from  knowledge 

of   shear   stress  (friction)  acting  on  the  fluid  adjacent  to   the 

wall.    Known  as  the  Reynolds-Colburn  Analogy,  it  is  described  by 

Rohsenow  in  the  following  manner  £2, p. 1853: 

"the  heat  transferred  to  a  surface  divided  by  the 

maximum   amount   which   could  be   transferred  in 

bringing   the   fluid  to  the   temperature   of  the 

surface    equals   the   momentum   transferred  in 
passing   over   a  surface  divided  by   the   maximum 

amount   of   momentum  which  can  be  transferred  in 
stopping  the  fluid  relative  to  the  surface." 

This  analogy,   exactly  true  in  parallel  flow  when  Pr  =  1   and 

reasonably   accurate  at  Pr  values  not  very  high  or  very   low, 


may   be  used  in  either  constant  heat  flux  or  wall  temperature 
conditions  and  may  be  expressed  mathematically  as: 

j  Colburn  -factor  =  StPr2/'-'  =  Cf/2  (9) 

where:   Stanton  # ,  St  =  h/pC_U  (10) 

C-f  =  Tw/<pU2/2)  =  -friction  coefficient      (11) 

Tw  =  apparent  shear  stress  acting  on  the  fluid 
at  the  wal 1 


Significant  -flow  experimentation  has  determined  that  -for  low 
turbulent  flow  conditions  over  a  smooth  plate  at  a  distance  x 
from  the  leading  edge: 

C.  =  .0592  Re  ~"2  (12) 

+  flat  plate  x 

Similarly,  for  turbulent  flow  conditions  in  a  cylindrical 
pipe  of  diameter  D: 

C,  =  .046  Ren~-2  (13) 

Tpipe  flow  u 

Substitution  of  these  relations  into  equation  (9)  allows  for 
the  determination  of  h: 


h,  ,  flat  plate  of  length  1: 


stPr2/3  _  c^/2 

hxPr2/3  =  .0592  (Re~'2)/2 

hx  =  .0296  pC  U(Rex)~-2  Pr~2/3 

=  .0296  pC  U(Rew)~-2  Pr~2/3  x  K 

K  x 
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substituting  -for  density,  the  equation  (6)  identity,  p  =  ^/v^..: 
hw  =  .0296  K  Ux  u.C  <Rew)~'2  Pr~2/3 

,-%  —  — —  i  .  ■  ■  p      yi 

X  vk   K 

recognising  equations  (4)  and  (5)  identities  Re„  =  Ux  and 

Pr  =  ^C  /K  : 


h.,  =  .0296  K  Re  '8  FT17 


integrating  over  length  1  to  determine  mean  h-. 


hl   "if1  -0296  ti  (UK/vk'*"B  Prl/3  dx 


Jo 


1   J,-,  X 


=  .0296  K  (U/t^.)'8  Pr1/3  f1  x  ,2  dx 
1       %  J0 

=  .0296  K  (U/v..)'8  Pr1/3  Cx"8!!1 

—      *•  o 

0.8   1 

h1   =  .037  ^(Re^'8  Pr1/3  (14) 

1 

A  similar  substitution  of  equation  (13)  -for  turbulent  -flow  in 

a  cylindrical  pipe  yields: 


hfv,  pipe  flow  o-f  diameter  D: 


hn  =  .023  K(Ren)'8  Pr1/3  (15) 


These  particular  relations  form  the  basis  of  many  heat 
transfer  correlations  commonly  used,  rigorously  applicable  in 
fully  developed,  steady  turbulent  flow  within  established  Re 
and  Pr  limits  where  temperature  difference,  AT,  across  the 
boundary  layer  is  not  extreme. 


2.  MODELING  CONVECTIVE  HEAT  TRANSFER  IN  RECIPROCATING  ENGINES 

In  reciprocating  engines,  the  significant  process  by 
which  termal  energy  is  lost  is  by  -Forced  convective  heat 
transfer  -from  the  gas  medium  across  cylinder  inner  surface 
boundaries.  Large  variations  in  the  magnitude  and  direction 
of  convective  heat  transfer  exist  from  engine  to  engine  and 
point  to  point  within  an  individual  engine.  The  conditions 
influencing  heat  transfer  Are  so  complex  and  locally 
different  through  an  engine  cycle  that  rigorous  mathematical 
interpretation  is  invalid.  Radical  variations  in 
instantaneous  flow  patterns  allow  gas  velocities  to  be  only 
vaguely  characterized.  Non-uniformity  of  physical  properties 
in  the  gas  medium  and  surface  materials  exists.  The 
combustion  process  is  particularly  difficult  to  model  and 
further  complicated  by  compression  ignition.  Two  of  the  most 
serious  obstacles  question  the  fundamental  assumption  of 
fully  developed,  steady  turbulent  flow  and  the  concept  of  a 
finite,  instantaneous  value  of  h. 

In  operating  engines,  the  entire  cycle  is  completed 
within  milliseconds,  making  accurate  measurement  of  local 
instantaneous  heat  flux,  q,  in  itself  quite  difficult. 
Turbulent  intensity  in  four  stroke,  CI  engine  cylinders  is 
produced  by  the  suction  velocity  of  the  induction  stroke,  the 
intensity  being  roughly  proportional  to  it.  To  optimize 
fuel-air  micro-mixing,  off-set  induction  ports  may  be 
incorporated  to  create  swirling  gas  motion.  In  the 
compression    stroke   within   a   simple   geometry    cylinder, 
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induction  turbulence  is  suppressed.  For  swirl  enhanced 
cylinders  with  combustion  chambers  cut  into  the  piston  crown, 
conservation  of  angular  momentum  in  the  cup  created  by 
viscous  drag  accelerates  swirl  to  significant  angular 
velocities  as  the  piston  reaches  TDC.  Radial  "squish" 
velocities  are  created  near  TDC  and  combine  with  swirl  to 
produce  a  complicated  toroidal  motion.  In  the  expansion 
stroke,  rapid  deaccel erati on  of  rotational  velocities  in 
direct  injection  cylinders  are  superimposed  by  i 1 1 -def ined , 
small  combustion  induced  velocities  as  hot  combustion  product 
gases  expand;  in  indirect  injection  cylinders,  prechambers 
deliberately  channel  and  exaggerate  the  charge  velocities  to 
*  blow  out  of  energy  cells  at  levels  up  to  500m/sec  [33.  In 
the  exhaust  stroke,  the  combustion  products  are  pushed  out  by 
the  piston.  The  accumulation  of  these  factors  indicate  that 
turbulent  boundary  layers  are  not  steady,  but  are 
periodically  created  and  destroyed  in  a  manner  beyond  our 
ability  to  model.  In  addition,  because  of  the  thermal 
capacity  of  the  gases,  there  is  a  thermodynamic  phase  lag 
markedly  apparent  in  reciprocating  engines;  AT  and  q  do  not 
pass  through  zero  at  the  same  instant  indicating  values  of 
h  =  ±co  and  highly  fluctuating  values  immediately  about  it. 

Recognizing  these  obstacles,  investigators  have 
concluded  that  within  obtainable  accuracy  it  is  still  best  to 
model  heat  transfer  as  a  quasi -steady  process  and  to 
correlate  cyclic  observations  in  terms  of  conventional  Nu , 
Re,    Pr   numbers.    Values   of   the   various   constants   and 
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exponents  Are    determined  in  individual  engine  tests;  a  wealth 

of   experimental   data   provides   guidance   -for    constant  (s) 

selection   for  operating  engines  or  design  studies  o-F  similar 

construction.   Hohenberg  states  C4,  pg.  2783]: 

"It  would  be  erroneous  to  look  -for  a  -formulation 
which  includes  all  possible  influences  on  the 
heat  transfer.  This  would  only  render  the 
calculation  more  complicated  rather  than  more 
accurate.  The  aim... must  be  to  find  out  and 
mathematically  describe  the  really  essential 
factors  influencing  heat  transfer  in  the  engine." 

If   the  time  for  change  in  bulk  stream  conditions  to   diffuse 

through    a   boundary   layer   is   small   relative   to    other 

significant   times  in  the  engine  cycle,   then  boundary  layers 

may   be   considered  quasi-steady  and  at  any  instant   be   that 

associated   with  existing  conditions.    Le  Feuvre  et   al   C53 

approximating  diffusing  time  as  &  '""i*  ana"  modeling  gas  motion 

in    solid    body   rotation,    successfully   performed    such 

calculations  in  a  direct  injection,  swirl  enhanced,  CI  engine 

operating   at  high  RPM  (2000  RPM) .    As  for  the  existence   of 

phase  lag,  Annand  C63  states  that  the  actual  quantity  of  heat 

transferred  during  this  period  is  minute;   the  error   arising 

from   the   use   of  a  finite  h  at  each  instant  of   the   piston 

stroke  is  negligible. 

The    majority    of   the   heat   transferred   within    a 

reciprocating    CI    engine   is   due   to    turbulent    forced 

convection.   On  the  basis  of  experimental  data  and  theoretical 

considerations,    this   relationship   is   expressed   by    the 

formula: 


Nu  =  constant (Re)n  Pr1  (16) 


Recognizing  that  Pr  x  0.7  -for  common  gas  mediums  due  to 
predominance  o-f  0--,  and  N-0 ,  substitution  of  the  ideal  gas  law 
identity,  p    =  F'/R'Tg,  into  an  expanded  equation  (16)  yields: 

h  =  constant  K^~nPnUnTg~nl cn_1  (17) 

Assuming  that  actual  gas  velocities  will  be  at  least 
approximately  proportional  to  mean  piston  speed,  a  -function 
o-f  Sp  is  used  with  the  multiplicity  constant  determined  by 
engine  design.  Tg  is  used  vice  the  average,  (Tg  +  Tw)/2,  as 
the  calculation  o-f  Tg  -from  measured  instantaneous  pressure 
and  gas  law  application  already  represents  an  average  -for  the 
cylinder  contents.  A  consideration  o-f  gas  disassociati  on  is 
neglected,  as  any  attempt  to  quantify  it  would  be  eclipsed  by 
the  previous  assumption  of  uniform  properties.  During  the 
combustion  process,  gas  product  composition  is  assumed  to 
equal  the  burned  composition.  Tables  and  Figures  of  thermal 
properties  (>».,  C  ,  Pr ,  3f)  for  individual  gases  and  CnH0  lean 
hydrocarbon  fuel  combustion  products  Are  provided  in  Appendix 
C  £6,73.  Cyclic  local  instantaneous  surface  temperature 
profiles,  from  which  actual  measurement  of  q  is  generally 
derived,  Are  demonstrated  by  Oguri  et  al  C83  to  show 
transient  magnitudes  typically  less  than  50°F  for  sooty 
surfaces  and  less  than  20°F  for  clean  surfaces.  With  the 
exception  of  exhaust  valve/port  entries  which  may  typically 
see  temperatures  greater  than  1000°F,  the  temperature  range 
exhibited  by  cylinder  sleeve,  head,  and  piston  crowns  are 
fairly  consistent  allowing  assumption  of  a  mean  Tw. 


3 .      WOSCHNI 'S   CORRELATION   FOR   MEAN   INSTANTANEOUS   HEAT 

TRANSFER  COEFFICIENT,  h„,-  IN  CI  ENGINES 

A  recognized  industry  standard  for  -four  stroke,  direct 
injection,  quiescent,  CI  engine  heat  transfer  calculation  is 
Woschni's  correlation,  stated  in  U.S.  units  as  follows  C9D: 

he  =  15.480  P*8CC1Sp+C2(VTgc/PcVc) (P-Po) D,8B~*2Tg~-53  (18) 
where  he  =  BTU/hr-f t2-0 R 

P  ,  V  ,  Tg   =  instantaneous  values  at  ignition  start 
Po  =  motoring  pressure 
Based   on   suppressed   swirl  engine  designs,   the   values   of 
constants  C<  and  C-,  s^rei 

1.  induction  stroke       C1  =6.18    C0  =  0   <18a) 

2.  compression  stroke     C.  =  2.28    C^  =  0   (18b) 

3.  expansion  stroke  (motoring)     C<  =  2.28 

C~  =  0 


(18c) 


4.  combustion-expansion  stroke     C..  =  2.28 

C^  =  5.9054  x  10~3  ft/sec-oR 


(18d) 


5.  exhaust  stroke  C.  =  6.13  C->  =  0  (18e) 
From  equation  (17),  a  n  exponent  of  0.8  was  determined  by 
Woschni  by  successively  changing  one  parameter  at  a  time  over 
maximum  possible  range.  A  consensus  opinion  of  all  recent 
published  correlations  is  that  0.7  =  n  =  0.8,  reflecting 
agreement  with  the  laws  of  similarity.  Woschni  developed  the 
motoring/firing  correlation  and  scavenging  correlation  on 
suppressed  swirl  design  cylinders  stating  C9, p. 30713s 
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"The  gas  motion  generated  during  the  induction 
phase  dies  away  with  time,  owing  to  internal 
-friction.  In  the  compression  stroke,  however, 
dependent  upon  the  piston  or  cylinder  cover 
shape,  an  additional  gas  motion  can  be 
induced  Csqui sh— swir 1 1  whose  intensity  is  also 
proportional  to  piston  speed.  In  this 
particular  case,  the  engine  involved  had 
direct  injection  and  straight  forward 
combustion  chamber  walls,  so  that  decay  due  to 
internal  friction  prevailed." 

A   Tg~"*uC'  exponent  is  due  to  Woschni's  assumption  of  air   and 

—       —   £>        —       —   75 
gas  product  >*.(Tg)  «  Tg   and  K(Tg)  «  Tg     at  n  exponent  =  0.8 

(see  equation  (17)).    To   approximate   ill-defined   velocity 

increase  due  to  combustion,  the  difference  P-Po  which,  unlike 

Tg— To,   sees   a   rapid  decay  upon  combustion   completion,   is 

used,   multiplied  by  instantaneous  cylinder   volume,   v" ,   and 

related  to  the  charge  mass  at  beginning  of  combustion, 

mc  =  PcVc/RTg_. 

It   is   at   this  point  that  Woschni   makes   a   critical 

simplifying  assumption.    Recognizing  the  significant   effect 

of   luminous   flame   radiation  heat  transfer  in   CI   engines, 

Woschni   arbitrarily   accounts   for  combustion   radiant   heat 

* 

transfer  by  affixing  the  value  C0  to  equate  actual  measured  Q 
integrated  over  the  stroke  with  Q  combustion  determined  by 
the  correlation.   Woschni  states  C9, p. 30753: 


"There  is  no  doubt  that  such  an  influence  occurs 
in  the  engine,  and  it  is  taken  into 
consideration  in  lumped  form  in  the  combustion 
term  as  defined  here,  but  it  seems  impossible  to 
the  author  to  separate  it  from  the  convective 
term.  The  only  question  that  remains  to  be 
considered,  therefore,  is  the  error  made  regard 
to  the  tendency. .. for  the  radiation." 
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In  essence,  for  application  to  the  New  Cycle  Engine 
expansion  cylinder  operating  at  significantly  lower  peak 
pressure  conditions,  laws  of  similarity  allow  significant 
extrapolation  of  Woschni's  correlation  for  forced  convection 
heat  transfer  only.  The  correlation  must  be  corrected  to 
account  for  the  fact  that  radiant  heat  transfer  is  not  a 
function  of  <P-Po)"8. 

4.        CALCULATION    QF   ORDER   OF   MAGNITUDE   REDUCTION   IN 
COMBUSTION  FORCED  CONVECTION  HEAT  FLUX 

The  representative  firing  cycle  P-Tg  versus  crank  angle 
plot  published  by  Woschni  for  natural  aspirating  operation  is 
shown  as  Figure  2—1  C93.  The  range  of  additional  principle 
operating  parameters  over  which  the  motoring/firing 
correlations  were  developed  well  representive  the  essential 
parameters  of  New  Cycle  Engine  <NCE)  expansion  cylinder 
operation  C93: 

Parameter     Noshni  (W1)  engine     NCE  Expansion  Cylinder 

RPM  900  max  818.3 

Sp        9  m/s  max  C29.6  f/sU       33.33  ft/sec 

Supercharge         3  max  1 . 0 

pressure 
ratio 

Tw        nominal  500° K  C900°R3  900° R 

L,B  unknown  1.22204  ft 

Bore   and  stroke  dimensions  were  not  provided,   however   from 

known   maximum  data  point  values  Sp  =  9  m/sec  and  Nppjv.  =  900, 

the  approximate  stroke  length  L  may  be  determined: 


■j.to 


L  SB    60    Sp  ss   (60)  (9)   as  0.3  m  CO.  98  -ft: 


">t 


-NRPM  <2) (900)  appro*. 

For  conventional  CI  engine  geometry,  bore  diameter  B  is 
approximately  equal  to  L;  there-fore  -for  effective  comparison 
of  he  proportional  to  B-  in  -forced  convection  heat 
transfer,  consideration  of  bore  diameters  of  similar 
magnitude  dimensions  may  be  neglected.  The  significant 
difference  between  combustion— expansi on  stroke  profiles  are: 

1.  An  order  of  15  magnitude  increase  in  peak 
combustion  pressure  in  Woschni's 
conventional  high  r  ,  CI  engine. 

2.  A   reduced  Tg  decay  rate,   thus  increased 

(Tg-Tw)  forced  convective  heat  transfer 
potential  in  the  New  Cycle  Engine 
expansion  cylinder. 

3.  A  delayed  fuel  injection/ignition  start 
in  the  New  Cycle  Engine  expansion 
cyl inder . 

Assuming  that  the  magnitude  of  mean  instantaneous  gas 
velocities  through  the  expansion  strokes  Are  roughly  equal, 
the  approximate  ratio  of  forced  convective  heat  transfer  loss 
in  the  respective  engines  may  by  obtained  by  determining  the 
mean  'ratio  of  instantaneous  heat  flux,  q,  over  a  90°  crank 
angle  domain  from  ignition  start.  A  selection  of  90°  is 
based  on: 

1.    Establishing  a  normalised  basis  in  which   to 
compare   the  order  of  magnitude  reduction  in 
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forced  convection  versus  radiant  heat 
transfer;  the  crank  angle  domain  over  which 
luminous  flame  radiation  is  significant  in 
medium  speed  CI  engines  less  than  1000  RPM 
is  nominally  90°  from  combustion  start  C8]. 
2.  At  90°  the  slope  of  pressure  and 
temperature  decay  transients  induced  by- 
combustion  in  both  engines  become  small. 

Setting  the  radiant  heat  transfer  constant  C^  =  0  in  equation 

(18d)  : 

qW/qNCE  °*  ^e(Tg"Tw)W/Ffe(Tg~Tw)NCE 


CP-8Sp*8/B-2Tg,53:CTg-TwDw 
qW/qNCE  « 


CP" 8Sp ' 8/B' 2Tg " 533  CTg-TwDNCE 


Neglecting   B  consideration  and  assuming  equal  magnitude   gas 
veloci  ties: 


CP"8/Tg"53]CTg-Tw:, 


qw         tF"-/rg---3Crg-TwJw 

«   (19) 

8—  ^f 
/-r  r-.   •  ««  I 


qNCE       CP'a/Tg---i:Tg-Tw:NCE 

where  Tw  =  900° R 

Table  2-1  summarizes  the  calculations  based  on  Figure  2-1  and 
Appendix  I  data  at  5°  increments  from  combustion  start.  The 
mean  heat  flux  reduction  ratio  may  then  be  determined  by 
numerical  integration  using  Simpson's  rule: 
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qw  =       _^_       f90    Cqw/qNCE:de        (20) 

qwrF  A90°     J,-. 

INLtl-forced    convection  u 


A5°         qW           qW           qW                     qw           qw           qw 
CI +4 +2 +...+2 +4 +1 


■UJ 


qNCE 


qNCE   qNCE   qNCE       qNCE   qNCE   qNCE 

80     85     90° 

A90° 


=   A5°/5  C328.673    =     6.087   approx . 
90° 


/ 
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TABLE    2-1 


Calculation  of  order  o-f  magnitude  reduction  in  combustion 
-forced  convection  q  -from  Woshni  '  s  test  engine  to  New  Cycle 
Engine       expansion       cylinder       operation       at       5°       crank  angle 

increments 


Crank  angle 


0 

5 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 

80 

85 

90 


Woshni 

engine  (M) 

p 
/c.2 

Tg 

P 
ati 

Tg 

•R 

40 

875 

38.71 

1575.0 

60 

1400 

58.07 

2520.0 

63 

1525 

60.97 

2745.0 

60 

1600 

58.07 

2880.0 

55 

1650 

53.24 

2970.0 

47 

1630 

45.49 

2934.0 

40 

1600 

38.71 

2880.0 

33 

1540 

31.94 

2772.0 

26 

1480 

25.16 

2664.0 

22.5 

1440 

21.18 

2592.0 

20 

1400 

19.36 

2520.0 

17 

1365 

16.45 

2430.0 

IS 

1330 

14.52 

2394.0 

12.5 

1305 

12.1 

2349.0 

11 

1280 

10.65 

2304.0 

10 

1255 

9.68 

2259.0 

9 

1230 

3.71 

2214.0 

B.5 

1215 

8.23 

2187.0 

8 

1200 

7.74 

2160.0 

NCE  Expansion  Cylinder 


p 

ati 

Tg 

•R 

4.15 

1896.6 

4.25 

2252.6 

4.09 

2571.2 

3.92 

2864.1 

3.76 

3136.0 

3.26 

3064.2 

2.89 

2998.7 

2.59 

2939.1 

2.34 

2884.8 

2.13 

2835.2 

1.96 

2790.0 

1.81 

2748.8 

1.69 

2711.3 

1.58 

2677.1 

1.49 

2646.0 

1.42 

2617.8 

1.35 

2592.3 

1.3 

2569.2 

1.25 

2548.4 

VqNCE 
equation  (19) 

4.46 

9.14 

9.29 

8.68 

7.94 

7.92 

7,68 

6.66 

6.2 

5.76 

5.65 

5.16 

4.92 

4.46 

4.17 

3.97 

3.75 

3.68 

3.59 


40 


80. 
(1600) 


rJP    60. 
^(1200) 


cu      40, 
£  (800) 

•  • 
QO 

£       20; 

(400) 


0 


360  400 

(TDC) 

0    I       degrees 


480 


FIGURE  2-1:  Woschni's  test  engine  firing  cycle;  natural 
aspiration  (_9J 
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5.      ESTIMATION  OF  MEAN  GAS  VELOCITY  MAGNITUDE, D 

5. 1    CONSTRUCTION  OF  CHARGE  IMPULSE  VELOCITY , Ui , TRANSIENT 

Acceptance  of  the  general  applicability  of  Woschni's 
correlation  to  the  New  Cycle  Engine  expansion  cylinder 
balances  on  the  assumption  o-f  comparable  gas  velocity- 
magnitudes.  While  geometric  flow  patterns  are  indeterminate 
in  preliminary  design,  the  magnitude  of  the  charge  impulse 
into  the  expansion  cylinder  leaving  the  regenerator  port  can 
be  ascertained.  Recognizing  the  inevitable  deaccel erati on  of 
this  impulse  in  a  simple  geometry  cylinder  on  the  expansion 
stroke,  a  conservative  estimation  of  the  mean  bulk  gas 
velocity  at  the  walls,  Uexp,  through  the  expansion  stroke  can 
be  made  with  minimum  error.  Substituting  Uexp  directly  into 
the  steady,  fully  developed  turbulent  flow  equations  (14)  and 
(15)  for  flat  plate  and  pipe  flow  should,  if  the  magnitudes 
of   cylinder  velocities  are  in  fact   similar,   be   reasonably 

* 

close  to  h  calculated  from  Woschni's  correlation.  To  reduce 
uncertainties,  the  relationships  will  be  first  compared  in 
motoring  conditions  from  TDC  to  45° ,  prior  to  ignition  start. 
Dent  et  al  C103  successfully  used  a  similar  method  in 
approximating  forced  convective  heat  transfer  loss  in  a  four 
stroke,  direct  injection,  high  swirl,  CI  engine  operating  at 
peak  TDC  swirl  ratios  equal  to  20.  Predictions  using 
Woschni's  motoring  correlations  failed  uniformly  low  compared 
to  measured,  spatially  averaged,  instantaneous  heat  flux. 
This  is  expected,  considering  the  central  assumption  of 
suppressed    swirl   inherent   in   application   of    Woschni's 
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correlation.    By  substituting  directly  into  the  flat  plate  h 
equation  (14),  measured  instantaneous  angular  gas  velocity  u> 
at  1   equal  to  bore  circumference  (2irr)  ,  remarkable  agreement 
was  obtained.   Agreement  was  not  as  good  at  -Firing  conditions 
due  to  the  transient  combustion  processes. 

From  ideal  gas  law,  at  a  given  angle  6  in  the  expansion 
stroke,  charge  mass,  m,  in  the  expansion  cylinder: 


Pv"      (2116.224  lbf /f t^-atm)PV 

m9  lbm  =  —   = 

RTg  (53.90)Tg 


=  39.262  PV  (21) 

Tg 

Taking  5°  crank  angle  as  the  unit  increment  of  time  at   818.3 

RPM;   the   charge  mass  rate,   entering  the  expansion  cylinder 

i  s: 


Am(5° ) 


m6  lbm/sec  ~    *"<5°>  = 


At  (5°  )     (NRp|v)/60)  360° 


Am(5°)  (22) 


.0010134  sec 


The  continuity  equation  for  impulse  flow  leaving  the 
regenerator  port  cross-section  area,  A  ,  equal  to  0.80  ft""" 
is: 


me  =  pUiAp  (23) 

Therefore,   by   rearranging  equation  (23),   impulse   velocity 
leaving  the  regenerator  port,  Ui  ,  can  be  determined: 


4? 


meRTg  me<53.90)Tg 

Ui-ft/sec  =  =  n 

PA      (2116.224  lb-f /ft*-atm)P(.080) 


31837  me  Tg 

(24) 


A  summary  o-f  the  calculations  using  instantaneous  P,  Tg ,  V 
values  -from  Appendix  A  at  5°  increments  is  provided  in  Table 
2—2.  Charge  mass  as  a  -function  o-f  expansion  stroke  crank 
angle  is  plotted  in  Figure  2—2,  indicating  that  the  main 
impulse  -from  the  compression  cylinder  via  the  regenerator 
exit  port  completes  10°  past  ignition  start. 

At  45°  combustion  start   "/-^total  =  2.1905x10""?  -  89.37. 


2. 4537 x 10  * 


At  55°  main  impulse  completion 


7-mtotal  =  2.28894x10^  =  93.3'/. 
2. 4537 x 10  * 

Ui   as  a  function  o-f  expansion  stroke  crank  angle  is   plotted 

in   Figure  2-3,   indicating  a  peak  Ui  equal  to  519.95   ft/sec 

near  -fuel  injection  start  (30°). 
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TABLE  2-2 

Calculation  o-f  charge  Hi  leaving  the  regenerator  port  into 
the  New  Cycle  Engine  expansion  cylinder  at  5°  crank  angle 
increments 


Crank  Angle  m 

8°  lbm 

eq. (21) 

TDC  0 

5  1 . 5550E-4 

10  7.7090E-4 

15  2.1309E-3 

20  4.5S21E-3 

25  8.3113E-3 

30  1  .  2936E-2 

35  1 . 7373E-2 

40  2.0491E-2 

45  2. 1905E-2 

50  2.2887E-2 

55  2.2889E-2 

60  2.2932E-2" 

65  2.2992E-2 

70  2.3056E-2 

75  2.3307E-2 

80  2.3512E-2 

85  2.3680E-2 

90  2.3820E-2 

95  2.3937E-2 

100  2.4035E-2 

105  2.4118E-2 

110  2.4188E-2 

115  2.4247E-2 

120  2.4298E-2 

125  2.4342E-2 

130  2.4379E-2 

135  2.4411E-2 

140  2.4438E-2 

145  2.4462E-2 

150  2.4480E-2 

155  2.4496E-2 

160  2.4509E-2 

165  2.4519E-2 

170  2.4527E-2 

175  2.4532E-2 

BDC  2.4537E-2 


Amc-o 

m 

Ui 

AlbmE-3 

1 bm/sec 

ft/sec 

eq. (22) 

eq. (24) 

0 

0 

0 

0 . 1 555 

0. 1527 

41.82 

0.6154 

0.6043 

133. 14 

1 .  360 

1.3275 

236.92 

2.451 

2.4069 

--•*-!.■£.  a  69 

3.730 

3.6618 

458. 13 

4.625 

4.5413 

519.95 

4.437 

4.357 

499. 15 

3.  118 

3.0614 

382.68 

1.414 

1 . 3882 

202.03 

0 .9816 

0.9641 

162.73 

0. 0028 

0. 00275 

0.55 

0 . 0426 

0.04183 

9.72 

0.0596 

0.05852 

15.54 

0.0641 

0.06294 

18.83 

0 . 25 1 0 

0.2465 

81.44 

0.2049 

0 . 20 1 2 

72.76 

0. 1687 

0. 1657 

65.  12 

0. 1397 

0. 1372 

58.  15 

0. 1167 

0. 1146 

52.03 

0 . 098 1 

0. 09633 

46.56 

0.0827 

0.08121 

41.54 

0.0702 

0.06893 

77   1  O 
•_'  /  .  1  ^_ 

0.0596 

0.05852 

33 .  0 

0 . 05 1 1 

0.05018 

29.84 

0.0436 

0.04281 

26.  13 

0.0372 

0.03653 

23.04 

0 . 0320 

0.03142 

20.41 

0 . 0270 

0.02651 

17.67 

0 . 0239 

0.02347 

15.99 

0 . 0 1 83 

0.01797 

12.48 

0.0158 

0.01552 

1 0 .  94 

0. 0130 

0.01276 

9.  12 

0 .  0 1 0 1 

0.00992 

7.  15 

0.0074 

0.00727 

5.28 

0 . 005 1 

0. 00501 

3 .  65 

0.0052 

0 . 005 1 1 

3.70 

.025 


.020 


.015- 


.010- 


005-- 


TDC 


0   !     °de9Pees 


FIGURE  2-2:  Expansion  cylinder  mass  charge,  m,  versus 
expansion  stroke  crank  angle 
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FIGURE  2-3:  Expansion  cylinder  chargeVi  versus  expansion 
stroke  crank  angle 
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VERIFICATION  OF  TURBULENT  FLOW 


To    veri-fy  turbulent   -flow   conditions   leaving    the 

regenerator  port,  calculated  Re  at  30°  using  1   equal  to  the 

regenerator   port  cross-section  minor  dimension   (depth)   of 
0.75  inch  is: 


Uil     Uil    «Uil     P 
Re300      =  =  = (25) 


vk     p./p  p.  RTg 


From  Appendix  C,  Table  C-l  ,  ,H.(Tg)  and  C  (T)  -for  air  is: 

>u  lbm/ft-sec   =   2. 163x 10~7 (Tg) * 645  (26) 

(CI) 

C   BTU/lbm-hr   =   0.224  +  0. 262x1 0~4  Tg       (27) 

(C2) 


Substituting   instantaneous  values  of  P,   Tg  from  Appendix   A 
and  >i(Tg)  identity  equation  (26)  into  equation  (25): 

P30o  =  6. 10171  atm 
Tg30o  =  2194. 29°R 

(519.95) (0.75/12) (2116.224  lbf /ft2-atm) (6. 10171) 


(2.163x10  ') (2194.29) * u^°    (53. 90) (2194. 29) 


=  1.1476x10-  turbulent 
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5.3    ESTIMATION  OF  EXPANSION  STROKE, Uexp 

The  calculation  o-f  the  mean  impulse  velocity,  Di  ,  -for 
the  purpose  o-f  Uexp  approximation,  is  correctly  restricted  to 
the  main  impulse  domain  TDC  to  55°  only,  neglecting  residual 
10"/.  charge  mass  regenerator  bleed.  Combustion  induced 
perturbations  during  crank  angles  45°-55°  may  be  neglected. 
Ui  may  then  be  calculated  using  numerical  integration  by 
trapezoidal  rule  and  is  equal  to: 

0i-ft/sec  =   -A-  f55°  Ui  de  (28) 

TDC-550    A55°JTDC 

(A5°/2)C  UiTDC  +2Ui50  +2Ui  1(-)0  .  .  . +2Ui  5(-)0+Ui550  3 

A55° 

=  (A5°/2) C5979.033   =  271.77  ft/sec 


A5 


"jo 


A  reasonable  assumption  must  be  made  as  to  the  expected 
value  o-f  Uexp  over  the  entire  stroke.  The  principle  -factors 
causing  charge  velocity  deaccelerati on  include  cross-section 
Area  expansion  exiting  the  regenerator  port,  internal 
-friction,  cylinder  volumetric  expansion,  and  angular  momentum 
deaccel  erati  on  o-f  any  existing  rotating  -flow  patterns. 
References  C5,  10,  11,  12,  133  unanimously  agree  that  in  swirl 
enhanced  engines,  peak  instantaneous  swirl  at  TDC  is  much 
greater  than  the  mean  swirl  valves  over  the  stroke. 
Le  Feuvre  et  al  C53  calculated  -for  a  4.5  inch  bore  direct 
injection,  CI  engine  operating  at  2000  RPM,  that  10"/.  velocity 
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decrease  occurred  each  stroke  due  to  viscous  dissipation 
alone.  A  conservative  estimate  of  Oexp  =  (Oi Tnc_550 ) /2  is 
therefore  assumed;  selection  of  this  value  should  reduce  any 
error  in  either  direction.  It  is  noted  that  Ui  ynr._onrj  "from 
Table  2-2  in  firing  conditions  equals  102.9  ft/sec. 
Combustion  induced  velocities  due  to  charge  gas  expansion  in 
spark  ignition  (SI)  engines  have  been  calculated  by  Annand 
C3J  to  be  on  the  order  of  33  ft/sec,  with  expected  velocities 
in  CI  engines  less  due  to  diffusion  flame  configuration. 
Thus,  let  Oexp  =  (Di ync-550 ) /2  =  135-9  ft/sec. 
6.     COMPARISON   QF  h"e  CALCULATED  FROM  WOSHNI'S  CORRELATION 

AND  THEORETICAL  MODELS  ' 
6.1    EXPANSION  STROKE;  MOTORING  CONDITION 

■  I       I,     ■»■ — ■     i     I,  i      i    

Using  the  above  first  approximation  of  Oexp,  estimates 
of  h  from  theoretical  equations  (14)  and  (15)  for  steady, 
fully  developed  flow  may  be  reached.  Expanding  pipe  flow 
equation  (15)  at  1   =  B  and  U  =  Oexp  yields: 

hn  =  .023(K/B)ReD"8Pr1/3  =  . 023 (K/B) C0expB/vk 3 ' BPr 1/3 

Substituting  the  identities  v..  =  /t/p  ,  K  =  C.^/Pr ,  and 

K  p 

P  =  P/RTg: 


hn  =  .  023C  u.CPUexp/^RTg3"8B  "  2Pr  2/3 


Substituting   the  C  (T)  and  ^(Tg)  identities   from   equations 
(27)  and  (28)  and  appropriate  conversions: 
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hD  =  . 023 C. 224+. 262x10  4Tg 3 C2.  163x 10  7 (Tg ) * 645 (3600sec/hr ) D* 
C (2116.224  lbf/ft2-atm  P  Uexp ) / (2. 163x 10~7Tg" 645R?g) 2 " 8* 
B  ~*Pr    *-' -* 

which  simplifies  at  known  R,  Uexp,  Pr  values  to: 

(135.9  P) "8 

hn  =  9. 194C2.24+2.62x 10  4Tg3 (29) 

BTU/hr— ft^-oR  B"*  Tg'*'1 


An   identical  expansion  of  the  -flat  plate  equation  (14)  at 
1   =  ttB  ,  expansion  cylinder  circumference,  yields: 

(135.9  P) ■ 1S 


hj,  -  1 1.  674 C 2.  24+2.  62x10  HTg  H —- —  (30) 

BTU/hr-ft^-oR  B*2  Tg"671 

A  summary  of  h  calculations  -for  the  charge  induction 
monitoring  angles,  TDC-450 ,  comparing  equations  (29)  and  (30) 
to  Woshni's  motoring  correlation  equation  (18c),  is  provided 
in  Table  2-3.  The  agreement  is  quite  good,  even  with 
acceptance  o-f  reasonable  error  in  Uexp  approximation. 


TABLE  2-3 

Comparison  of  R  calculated  by  Woshni's  correlation  and 
theoretical  models  during  charge  induction  into  the  New  Cycle 
Engine  expansion  cylinder 

Crank  Angle  h"e         hn        hj^ 

9°  eq.<13c)    eq.  <29)    eq.(30) 

5  16.994  15.221  19.477 

10  20.191  18.143  23.215 

15  24.006  21.593  27.631 

20  28.209  25.396  32.496 

25  32.029  28.839  36.902 

30  34.236  30.793  39.403 

35  33.S88  30.414  38.918 

40  31.168  27.897  35.697 

45  27.164  24.251  31.037 

h  in  BTU/hr— ft2-°R 
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6.2    EXHAUST  STROKE 

Attempts  to  directly  estimate  the  mean  bulk  gas 
velocity  at  the  walls  during  the  exhaust  stroke,  Uex h ,  in 
light  of  preliminary  design  restrictions  and  post-combustion 
ambiguities,  Are    not  valid. 

An  indirect  method  allowing  h  calculation  is  to  correct 
equation  (29)  by  a  factor  equivalent  to  the 
scavenging/motoring  Uexh  correction  observed  by  Woschni  -for 
suppressed  swirl  geometry,  direct  injection,  CI  engine 
cylinders.   From  equations  18a  and  13c: 

Uexh  correction  factor  =  (6.18  Sp/2.28  Sp ) ■ 8  =  2.22 

When  the  regenerator  exit  valve  3  opens,  the 
significant  flow  pattern  will  be  parallel  to  the  bore  axis, 
up  the  sleeve,  and  through  the  regenerator  port;  flow 
direction  is  therefore  analogous  to  pipe  flow.  From  Appendix 
C,  the  identities  ^(Tg),  C  (Tg),  Pr (J) ,  and  ^(Tg)  for  lean 
C  H~  hydrocarbon  fuel  combustion  products  at  P  =  1  atm, 
$  =  0.75  may  be  calculated: 


1.47x10  7(Tg)'7    1.47x10  7(Tg)*7 


>l  lbm/ft-sec 


1  +  .027$         1  +  .027 (.75) 

=  1.441x 10~7(Tg) -7  (31) 

(C3) 

C       BTU/lbm-°R    =    0.247    +    0. 243x1 0~4    Tg  (32) 

P  (C5) 

Pr<3T)     ■    0.05    +    4.2(2-1)     -    6.7Q-1)2  (33) 

(C4) 

3%(Tg)     z    Appendix    C,    Figure    C-l  (34) 


Applying  the  Dexh  correction  factor  =  2.22  and  the  identities 
^(Tg),   C  (Tg)   de-fined   above   to  pipe   -flow   equation   (29) 
yi  elds: 

hn  =  2. 22 (.023) (.247+.243X 10~4Tg) ( 1 . 441 x 10~7Tg " 7) (3600sec/hr ) * 


[(2116.224  lb-f/-ft2-atm  P  Uexp  )  /  ( 1 .  441  x  10  7Tg  *  7RTg  )  1 '  Q* 


">  /T 


g-.  *.pr  2/_. 

which  simplifies  at  known  R   ,  Dexp  values  to: 

135.9  P 

hD  =  14.835C2.47  +  2.43x10  4Tg 3 Pr  2/3 

BTU/hr— ft2-°R  B"2Tg'66 

(35) 

From  Appendix  A;  P  =  1.05  atm  and  2344. 53° R  i  Tg  i  2439. 06° R 

TDC  BDC 

through   the   exhaust  stroke.    A  comparison  o-f  h   calculated 

■from   equation  (35)  and  Woschni's  exhaust  stroke   correlation 

equation  (18e)  follows: 

crank  angle      Tg     <f(Tg)      Pr<3*)      he         hD 

°R(°K)   eq. (34)   eq.(33)   eq.(18e)   eq.(35) 

BDC       2439.06    1.275     0.698     17.587    17.062 
(1355) 

TDC       2344.58    1.280     0.701     17.959    17.331 
(1303) 

h  in  BTU/hr-f t2-°R 

The   above   series   of  calculations  suggest  that  to   a   first 

approximation,    application   of   Woschni's   correlation    in 

determining   forced  convective  heat  transfer  loss  for  the  New 

Cycle  Engine  expansion  is  valid. 


CHAPTER  THREE 

CALCULATION   QF   UNCORRECTED  HEAT  TRANSFER 
LOSS/HEAT  RELEASE  PER  CYCLE,  ft 

--■■-■-  ■     ■  i     i   '  i  i 

Table  3.1  summarizes  the  calculation  of  total 
instantaneous  heat  transfer  rate,  Q,  through  a  complete  360° 
-firing  cycle  o-f  the  New  Cycle  Engine  expansion  cylinder  using 
Woschni's  correlation  -for  h  .  For  calculation  o-f  h  in 
combustion,  motoring  P  is  generated  graphically  in  Appendix 
A,  Figure  A-l: 

G>BTU/sec  =  h"eAHT  (Tg-Tw)  ( 1  hr/3600  sec)   (36) 

Due  to  the  establishment  in  preliminary  design  o-f  a  simple 
cylinder  geometry,  A,,y  as  a  -function  o-f  crank  angle  may  be 
determined  directly  -from  instantaneous  V  values;  de-fining 
instantaneous  gas  volume  height  as  L^  and  neglecting  cylinder 
sleeve  clearance: 


E*2  Le 


AHT  =  2nB^    +    irBLe  =  2irB""  +  TrDC4V/irB*-3 


which  at  B=  1.22204  ft: 

AHJ   ^  =  2.346  +  3.273CV3  (37) 

ft2 


The  distribution  of  h  and  Q  as  a  function  of  crank  angle  is 
plotted  in  Figure  3-1.  It  is  immediately  clear  studying  the 
h  distribution  over  the  combustion  stroke  55°  to  BDC ,  that 
Woschni's  correlation  fails  to  correctly  extrapolate  luminous 


flame  radiant  heat  transfer  loss  to  low  r  operating 
conditions.  No  increase  in  h  occurs  due  to  the  predominant 
influence  of  increased  temperature:  associated  gas  layer  p 
decrease  and  )x.  increase  overcame  the  effect  of  K  increase, 
thus  lowering  h  ,  modeled  by  laws  of  similarity  solely  on 
forced  convective  heat  transfer.  The  low  r  and  delayed  fuel 
injection  operating  conditions  effectively  eliminate  any 
remarkable  pressure  rise  in  the  cylinder,  already  at  45° 
down-stroke  at  ignition  start.  Recalling  Woschni's  coupling 
of  radiant  heat  transfer  in  combustion  to  C/?(P-P  )  ,  the 
desired  "lumped  form"  term  effectively  disappears. 

Uncorrected  total  heat  transfer  loss  per  cycle,   Qj,  is 
determined   by  numerical  integration  under  the  0.  curve  using 
Simpson's  rule: 

QT  =  .00 10184  [DpDC+  4Q1850+  2Q19f)0  +  ... 

BTU/cycle       3 

+2Q170O+4Q175o+QBDC3      (38) 

=  . 0010184C7141.6133  =  2.42434  BTU/cycle 


From  Table  2-2,  total  charge  mass  per  cycle  equals  .024537 
lbm.  At  known  specific  heat  release  per  cycle,  A  heat,  equal 
to  500  BTU/lbm  charge,  the  ratio  of  total  heat  transfer 
loss/heat  release  per  cycle,  TU- ,  equals: 


UT  =  2.42434  BTU 


uncorrected    (500BTU/1 bm)  (.  024537  lbm/cycle) 


=  19.761"/. 


JO 


The  inability  o-f  Woshni's  correlation  to  account  -for  radiant 
heat  transfer  in  low  r  operation  is  clearly  demonstrated  by 
repeating  the  numerical  integration  with  the  combustion  term 
in  equation  (18d)  eliminated  (C->  =  0)  : 

Simpson's  Multiplier  =  7038.088 

QT  =  2.3892  BTU/cycle 

1lT  =  19.4747. 

uncorrected 

Considering  the  accuracy  of  the  correlation  in   general,   -for 

preliminary   design   prediction   purposes  setting  C-,  =   0   is 

considered  valid. 

To   correct  -for  luminous  flame  radiation  heat   transfer 


loss,  a  value  of  ^m^mpf   radi  ation  must  be  determined  which, 


in  comparison  of  ^w^NCE  -fQrcec'  convection  calculated  in 
Chapter  Two,  will  allow  first  approximation  of  the  required 
Uj   correction  factor.    As  there  exists  no  theoretical  basis 

for   coupling   magnitude  of  luminous  flame  radiation   to 

p 
(P— P  ) "  ,  an  attempt  to  determine  the  desired  ratio  from 

g 

(P-P  ) "  profiles  must  be  disregarded  in  favor  of  calculation 
founded  on  known  principles  of  radiation  heat  transfer  in 
diffusion  flame  combustion. 


vJ  / 


BTU 

AHT 
-ft-1" 

hr— ft2-°R 

TABLE  5-1 

Calculation   o-f   New  Cycle  Engine  expansion   cylinder   firing 

cycle   Q   using   Woschni 's   correlation  at   5°   crank   angle 
increments 

crank  angle           hB        AWT         Tg~Tw  Q 

6°  BTU 

sec 

eg  < 1 8 )    eg ( 37 )  Tw=900° R  eg ( 36 ) 

charge  induction eguati  on  ( 16c) 

5                16.944     2.357  1249.04  13.856 

10               20.191     2.390  1259.25  16.878 

16                24.006     "2.445  1295.73  21.128 

20               23.209     2.522  1325.26  26.190 

25               32.029     2.618  1329.71  30.972 

30               34.236     2.734  1294.29  33.652 

35               33.888     2.867  1216.61  32.834 

40               31.168     3.017  1110.57  29.009 

45               27.164     3.181  996.6  23.924 

combustion-expansi  on eguati  on  <  18d  ) 

50               26.086     3.358  1351.6  32.888 

55               23.969     3.548  1671.21  39.445 

60               22.066     3.741  1964.13  45.038 

65               20.853     3.944  2236.02  51.083 

70               18.710     4.152  2164.15  46.700 
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TABLE  3-1  (CONT'D) 


crank  angle 


y 


75 

80 

85 

90 

95 

100 

105 

110 

115 

120 

125 

130 

135 

140 

145 

150 

155 

160 

165 

170 

175 

BDC 


*e 

AHT 

Tg-Tw 

Q 

BTU 

°R 

BTU 

hr-f t2-°R 

sec 

1 7 . 084 

4 .  362 

2098.73 

43.444 

15.742 

4.573 

2039. 11 

40.775 

14.632 

4.783 

1984.75 

38.584 

13.715 

4.990 

1935. 18 

36.7S9 

12.841 

5.  192 

1889.98 

35. 002 

12.099 

5.388 

1848. 18 

33.478 

11.457 

5.577 

1811.25 

32. 148 

10.865 

5.757 

1777.09 

30.877 

10.438 

5.927 

1746.03 

30. 006 

1 0 .015 

6.087 

1717.83 

29.089 

9.652 

6.236 

1692.29 

28 . 36 1 

9.329 

1669.22 

27.567 

9.052 

6.498 

1648.44 

26.934 

8 .  830 

6.611 

1629.82 

26.428 

8.632 

6.711 

1613.21 

25.959 

8.467 

6.797 

1598.51 

25.554 

8.334 

6.871 

1585.62 

25.221 

8.229 

6 .  930 

1574.46 

24.941 

8.  151 

6.973 

1564.96 

24.725 

8.086 

7.011 

1557.07 

24.517 

8.044 

7 .  03 1 

1550.74 

24.363 

8 .  028 

7.037 

1545.46 

24.260 
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TABLE    3-1      (CONT'D) 


crank  angle 

Re 

AHT 

Tg-Tw 

Q 

0 

e 

BTU 

-ft^ 

°R 

BTU 

<— i 

sec 

hr~ft^-°R 

exhaust 

equation  ( 18e) 

BDC 

17.587 

7.037 

1539.06 

52 . 909 

185 

17.604 

7.031 

1534.63 

52.763 

190 

17.619 

7.011 

1530. 70 

cr  <->   cr  <-»  t 
J^.i  <_f ji Oi 

195 

1 7 . 634 

6.973 

1526.81 

52. 137 

200 

17.649 

6.930 

1522.95 

51.741 

205 

17.664 

6.870 

1519. 13 

51.208 

210 

17.678 

6.797 

1515.36 

50.578 

215 

1 7 . 693 

6.710 

1511.61 

49.394 

220 

17.707 

6.611 

1507.99 

49.035 

17.721 

6.498 

1 504 . 4 1 

43. 121 

x-  •_•  y 

17.736 

6.373 

1500.90 

47. 125 

235 

17.748 

6  •  j^-J'O 

1497.47 

46.037 

240 

17.761 

6.087 

1497.47 

44 . 870 

245 

17.774 

5.927 

1 490 . 89 

43.623 

250 

17.786 

5.756 

1487.75 

42.308 

255 

17.798 

5.576 

1484.71 

40 . 929 

260 

17.810 

5.388 

1481.71 

39.493 

265 

17.821 

5.  192 

1478.95 

38.012 

270 

17.832 

4.990 

1476.24 

36.488 

275 

17.842 

4.783 

1473.65 

34.933 

280 

17.852 

4.573 

1471. 18 
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TABLE    3-1      (CONT'D) 


crank  angle 

Re 

AHT 

Tg-Tw 

Q 

0 

e 

BTU 

■ft2 

OR 

BTU 

hr— ft2-°R 

sec 

285 

17.861 

4.362 

1468.82 

31.786 

290 

17.870 

4.  152 

1466.58 

30.226 

295 

17.879 

3.944 

1464.46 

28.685 

300 

17.887 

3.741 

1462.44 

27. 198 

305 

17.895 

3.545 

1460.54 

<->c?      —t-T—i 

310 

17.902 

3.356 

1458.73 

24.346 

315 

17.909 

3.  181 

1457.03 

23 . 057 

320 

17.915 

3.017 

1455.41 

21.851 

17.921 

2.867 

1453.88 

20.750 

330 

17.927 

2.734 

1452.41 

19.774 

335 

17.933 

2.618 

1451.01 

18.923 

340 

1 7 . 938 

1449.67 

18.217 

345 

17.944 

2.445 

1448.36 

17.651 

350 

17.949 

2.390 

1447.09 

17.244 

355 

17.954 

2.357 

1445.83 

16.996 

TDC 

17.959 

2.346 

1444.58 

1 6 . 906 

61 
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FIGURE  3-1 : 


Expansion  cylinder  firing  cycle  he  and  Q  versus 
crank  angle;  Woschni's  coorelation 
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CHAPTER  FOUR 

RADIANT  HEAT  TRANSFER  IN  THE  NEW  CYCLE  ENGINE 

EXPANSION  CYLINDER 

1.      THEORETICAL  BACKGROUND 

When  a  temperature  difference  exists  between  two 
bodies,  part  of  the  thermal  energy  is  transferred  as 
electromagnetic  waves  of  radiation  at  varying  wavelengths.  A, 
characteristic  of  the  body  temperature  and  identity.  In  the 
electromagnetic  spectrum,  thermal  energy  is  propagated 
through  .01  >im  £  A  ^  100  }un.  Bodies  above  absolute  zero 
temperature  continuously  emit  energy  which,  when  striking 
another  body,  may  be  partly  absorbed,  fo,  reflected,  ^,  or 
transmitted  ,  T,  where: 

jo  +  $  +  X  =  1  (39) 

Thermal  radiation  reflected  from  and  transmitted  through  a 
body  is  again  absorbed,  reflected  by,  or  transmitted  through 
any  additional  bodies  ultimately  establishing  an  equilibrium 
state.  An  ideal  body  that  absorbs  all  radiation  is  known  as 
a  black  body;  if  no  radiation  is  transmitted,  as  in  the  case 
for  most  industrial  liquids  and  materials,  the  body  is  known 
as  opaque.  Combustion  chambers  covered  with  a  thin  film  of 
oil  contaminated  with  soot  particles  will  absorb  a 
significant  percentage  of  incident  radiation,  T  =  0. 

Thermal  energy  emitted  from  a  body  unit  surface  per 
unit   time   at   a  particular   wavelength,   the   monochromatic 


emissive      power,       E^  ,       can       be    determined    as      a      -function      of 
temperature    T(°K)     by    Planck's    Law: 

E,  =    2Trccin2x~5Cccin/ecbxT-l]~1  (40) 

A    /       3        P        P 
erg/sec-cm 


where:      c  =     speed  o-f  light  2. 997902x 1010  cm/sec 

c  =   Planck's  constant  6.623777x10     erg-sec 

cfa  =    Boltzmann's  constant  1.38026x10""    erg/°K 

X  =     wavelength  in  cm 

T  =     temperature  in  °K 
E.   will   increase  with  increasing  A  to  a  maximum  value   then 

decreases,   the  distribution   shifting   toward   the   shorter 

wavelengths   at  higher  body   temperatures.   The   relationship 

between   T(°K)   and   Am,„   is  determined   by   Wein's   Law   of 

max  7 

Di  splacement : 

Xmau  T(°K)  =  2898  cm-°K 
max 

The  total  radiant  emissive  power,  E,  may  therefore  be 
calculated  at  a  given  body  temperature  by  integrating  the 
area  under  the  curve: 

E  =  f03  E^  dA  (41) 

Jo  x 

For  a  given  black  body,  substitution  of  equation  (40)  into 
equation  (41)  yields  the  Stephan  Boltzmann  Law: 

EBlack  Bodx   =  croT  (°R>  <42) 

BTU/hr-f t2 

where:  <rQ    =    Stephan  Boltzmann  constant  1 .  713x  10~9BTU/hr-f  t2°R4 
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Real  bodies,  however,  are  not  absolutely  black.  Those 
bodies  of  a  continuous  radiating  electromagnetic  spectrum 
where  the  E^  curve  is  proportional  to  that  o-f  an  absolute 
black  body  at  each  A  are  known  as  grey  bodies.  The  ratio  of 
emissive  power  o-f  a  grey  body  to  an  absolute  black  body  is 
known  as  the  coe-f  -f  ici  ent  of  emissivity,  £: 


E  grey  body  (T)    =  E  (43) 

E  black  body  (T) 


From  an  analysis  o-f  eguilibrium  emittance  and  absorptance 
between  a  grey  body  and  its  surroundings  at  the  trivial  case 
o-f  each  equal  to  the  same  temperature,  it  can  be  demonstrated 
that  £  is  equal  to  to  (Kirchoff 's  Law).  This  important 
phenomenon,  strictly  correct  only  within  the  de-fined 
limitations,  is  -fortunately  observed  to  be  approximately 
correct  in  many  real  conditions  where  temperature  differences 
are  not  severe.  Additionally,  many  real  industrial  materials 
radiate  near  grey  at  very  high  values  of  fc;  instantaneous 
radiant  heat  transfer  rate,  Qp  ,  between  two  near-grey  bodies 
1  and  2  of  different  temperatures  at  a  distance  apart  may 
therefore  be  calculated  from  equation  (44): 

QP1  -»  =8.  ^*nF.    oAUT(T14-To4)  (1  hr/3600  sec) 

Btu/sec 

(44) 


where:  £.,-,=   1  +  1-1 

1     "7 

F1  r-,    =  geometric  view  factor  between  bodies  1  and  2 


Tables  providing  -formulas  for  £  and  F  calculation 
between  multiple  bodies  at  various  common  geometries  are 
published  in  standard  heat  transfer  texts.  For  gas  or 
incandescent  solid  particles  radiating  as  a  defined  unit  body 
within  an  encompassing  simple  geometry  combustion  cylinder,  F 
may  be  assumed  equal  to  unity. 
2.  NON-LUMINOUS  GAS  RADIATION 
2. 1    CHARACTERISTIC  HEAT  TRANSFER  EQUATIONS 

Unlike  solid  bodies  which  radiate  across  continuous 
wavelength  spectrums,  gases  emit  and  absorb  energy  at  direct 
wavelengths  only  in  line  spectra.  For  simple  monatomic  and 
diatomic  gases  such  as  0^,N^»,  and  H0 ,  the  line  spectra  &rs 
negligible  and  may  be  considered  transparent  <T  =  1).  For 
complex  gases  such  as  generated  in  exhaust  products,  however, 
radiant  capacity  may  be  significant.  C0o  and  H^O  vapor 
predominate  as  radiant  emitting  gases  in  combustion  products; 
the  significant  bands  of  radiation  are  C14D: 

CO,-,  H^O 


Band 

1 

2.36 

3 .  02    /un 

Band 

<-> 

mm 

4.01 

4.80    hjh 

Band 

~9 

12.5 

-    16.5       >un 

*■>     o> 


1.24    -       3.27    /un 
4.8      -       8.5      >im 
12  -    25  /un 


In  radiating  gas,  both  emittance  and  absorptance  exist 
throughout  the  gas  volume.  Penetrating  radiant  energy  will  be 
reduced  by  absorption  proportional  to  the  number  of  molecules 
in  the  ray  path  as  a  direct  function  of  individual  gas 
partial  pressure,  temperature,  and  path  length,  X.  This 
intensity  attenuation  is  expressed  by  Beers  Law: 
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to^  ■  1-exp  <kAX)  (45) 

where:  k.  =  monochromatic  absorption  coefficient  per 

unit  length 

X  =  radiant  path  length 

A   gas   volume   is,    of   course,   not   grey.    For   ease   of 

calculation  an  equivalent  gas  emissivity,   £  ,   is  de-Fined  in 

accordance  with  equation  (46): 

£g  =  1  -  exp(kglR)  (46) 

where:  k   =  apparent  gas  absorption   coefficient   per   unit 
1 ength 

1  p,  =  effective  mean  beam  length 

For  a  defined  gas  volume  of  V  and  boundary  surface  A,,-r- ,  1R 
may  be  determined  from  equation  (47): 

1R    =  3.4V  (47) 

ft     AHT 

Qp.   between  the  gas  and  surrounding  wall  of  emissivity  £    may 
then  be  determined  from  equation  (48): 

QR         =  £  £wcrQAHTCTg4-Tw43  (1  hr/3600  sec)    (48) 
BTU/sec 

For  gas  volumes  containing  more  than  one  radiating  gas,  total 
£  is  calculated  by  addition  of  individual  gas  emissivities 
and  an  additional  negative  correction  factor,  -££  ,  for 
radiant  band  coincidence  which  reduces  the  sum.  For  internal 
combustion   engine   exhaust   gas   products   from   hydrocarbon 
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based   fuels,   the   resultant   total  gas   emissivity   may   be 
considered  equal  to: 


£   =  G   ■  +  G    -&£„  (49) 

9     CCU    H^O    9 


Again,   significant   data   base   has  allowed   development   of 

graphs  for  E    ,   E    ,   and  -AG,,  as  a  function  of  P ,   Tg  and 
CDn    H/-iU         y 

lp,  and  are  published  in  standard  heat  transfer  texts. 
Empirical  formulas  based  on  experimental  data  for  calculation 
of  emissive  power  of  C0-,  and  f-UQ  vapor  useful  for  a  first 
approximation  in  U.S.  units  are    C14D: 

Ern   =  1.7125xlO_5(Prn  1R) 1/3<Tg/100)3"5      (50) 


'CO   "   '  '  Xx"-"*  w   srC0  XR 
BTU/sec-ft2 


2 


"O  ra      >Di   ■  6  /  -i-_  /  j  i-w-i  \  •-'  •  U 


EH  Q  =  1.6738x10  °  Pw  n,alp,u(Tg/100)-',v      (51) 
2 
BTU/sec-ft2 


2.2    PREDICTION  OP  GAS  RADIATION  HEAT  TRANSFER 

References  C3,  6,  8,  9,  10,  14,  15,  16]  state  that  gas 
radiation  is  insignificant  in  calculating  high  r  CI  engine 
heat  transfer  loss.  To  ascertain  if  separate  consideration 
of  gas  radiation  in  the  low  r  New  Cycle  Engine  expansion 
cylinder  may  be  required,  it  is  sufficient  to  estimate  Qp 
magnitude  at  maximum  1R  and  Tg  in  the  exhaust  stroke  and 
compare  to  Q  calculated  by  Woschni 's  correlation.  For  a 
typical  diesel  fuel  of  C/H  mass  ratio  equal  to  6  operating  at 
lean  $,  distribution  of  combustion  product  gas  includes 
approximately    10%   CO-,   and   107.   H-,0   water   vapor    C141. 
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Calculation  of  QR  using  equation  <48)  and  (49)  at  Appendix  A 
parameter  values  at  exhaust  stroke  BDC  follows: 

Q  Woschni  =  52.909  BTU/sec 

V  =   1.43333  -ft3 

AHT         =   7.037  -ft2 

Tg         =  2439. 06° R 

P  =  1.05  atm:   F'H  Q  =  .105  atm      PCQ   =  .105  atm 

1R  =   3.4(1.43333)    =   .6925  -Ft 
7 .  037 

QCQ   as  (1. 7125x 10"5) C (. 105) (.6925) 3  * /3 ,2439.  06, 3' 5 (7. 037) 
2  i  00 

=  3.604  BTU/sec 

QH  0  x    (1.6738xl0~5)  (. 105) '8(.6925) ' 6 (2439. 06)  Z"  ° (7. 037) 
2  1 00 

=  0.226  BTU/sec 


QR    S£  3.830  BTU/sec 
gas 


In  the  above  calculations,  corrections  -for  A£Q?  £w,  and 
emissive  power  o-f  the  cylinder  surface  proportional  to  tr  Tw  , 
all  which  would  reduce  the  total,  have  been  ignored.  The 
estimated  value  equal  to  3.83  BTU/sec  is  well  absorbed  into 
52.909  BTU/sec.  Noting  that  for  all  other  crank  angles  in 
the  exhaust  stroke  both  Tg  an  1 p  decrease,  plus  in  the 
combustion-expansion  stroke  luminous  flame  radiation  exists, 
it  may  be  assumed  that  separate  consideration  of  gas 
radiation  in  the  New  Cycle  Engine  expansion  cylinder  heat 
transfer  calculation  is  not  required. 
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LUMINOUS  FLAME  RADIATION 


3. 1    FUNDAMENTAL  CAUSE 


The  phenomenon  of  -flame  radiation  is  quite  complex, 
dependant  intimately  on  the  physical  and  chemical  properties 
of  combustion,  pressure,  temperature,  and  a  myriad  of  engine 
and  -fuel  parameter  variables.  Two  basic  -flame  types  are 
produced  with  respect  to  radiant  capacity:  non-luminous  -flame 
characteristic  o-f  SI  engines,  and  luminous  -flame 
characteristic  o-f  CI  engine  di-f-fusion  -flames.  In  SI  engines, 
pre— mixed  -fuel  and  air  are  adi  abatical  1  y  heated  in 
compression  then  ignited  by  a  spark  resulting  in  combustion 
by  turbulent  propagation  o-f  -flame  through  the  charge  from  the 
point  of  ignition.  The  mixture  burned  at  the  beginning  of 
flame  travel  expands,  then  recompresses ,  reaching 
temperatures  hundreds  of  degrees  higher  than  that  which  is 
burned  last.  This  recompression  of  combustion  product  gases, 
particularly  CO/-,,  yields  visible  yellow  radiation.  No 
intermediate  solid  particles  are  formed,  however,  and 
apparent   emissivities   of   combustion   product   gases   as   a 

function   of  near  transparent  A  bands  is  small.    The  driving 

4 
potential   for   radiant  heat  transfer,   proportional   to   T  , 

diminishes  more  rapidly  then  forced  convection  heat   transfer 

proportional  to  T;   the  contribution  of  instantaneous  radiant 

heat   transfer,   QR,   integrated   over   the  engine   cycle   is 

i  nsigni  f icant . 

In    CI   engines,    the   fuel   spray    injecting    into 

adiabatical ly   compressed   air  near  piston  TDC  is  ignited   by 
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spontaneous  explosion  of  a  micro-mixture  of  fuel  and 
entrained  air,  thereafter  burning  diffusively.  There  Ar& 
many  sources  of  ignition  and  flame  tends  to  conform  to  the 
fuel  spray  shape  until  dispersed  by  gas  motion.  In  the 
chemical  combustion  of  hydrocarbon  based  fuels,  cracking  to 
simpler  structures  and  polymerization  occurs.  When 
insufficient  local  Q-.  exists  to  feed  the  rapidly  increasing 
reactivity  of  the  mixture,  thermal  decomposition  of  C-H  and 
C-C  bonds  occurs  resulting  in  incandescent  carbon  soot 
formation  radiating  in  a  continuous  spectrum.  Alcock  et  al 
C11H,  in  a  classic  high  speed  photographies  study  of  a  Meurer 
direct  injection,  medium  speed  CI  engine,  clearly  observes  a 
luminous  flame  front  filling  the  chamber  within  20°  of 
ignition  start  (10°  BTDC) ,  glowing  white  hot  at  peak  flame 
temperatures  equal  to  4500°R,  then  cooling  by  expansion  to 
dull  red.  Unlike  SI  engines  where  combustion  completes  early 
in  the  expansion  stroke,  the  CI  combustion  process  extends 
late  into  the  expansion  stroke. 

This  difference  between  local  0o  concentration  and 
cracking  rate  is  the  principle  cause  of  soot  formation  and 
luminous  flame  radiation.  The  magnitude  of  associated  heat 
transfer  cannot  be  ignored,  and  has  been  measured  by  Annand 
C33,  Sitkei  £141,  Dent  et  al  C103,  Flynn  et  ai  C153,  and 
Oguri  et  al  C83  as  typically  207.  to  30"/.  of  total  heat 
transfer  loss,  Qj  ,  in  non-swirl  enhanced  four  stroke,  direct 
injection,  high  r  CI  engines  under  varying  design  and 
operational  conditions.    Oguri  et  al  131    demonstrate  this  in 
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Figure  4-1,   showing  GU/Q-r  between  two  such  different  size  CI 

engines   operating  at  varying  RPM ,   load,   and  cylinder   wall 

soot   contamination:    Engine   A  is  a  5.9  inch   bore   natural 

aspirating  engine;   Engine  B  is  a  12.2  inch  bore  supercharged 

engine.   In   very  large  diesel  engines  the  contribution  to  GU 

of   radiant   heat  transfer  is   even   greater;   J.   Butler   of 

Dexford  Co.,   United  Kingdom,  estimated  the  radiant  component 

of   Qj   as  75"/.  for  a  40  inch  bore  diesel  C83.    In  any   case, 

measured   peak   flame   emi ssi vi ties,    £,,   approaches   unity 

(1.0),   and   while  reaction  zone  temperature  is  difficult   to 

determine,   a   consensus  opinion  shared  is  that  peak   radiant 

temperature  is  much  higher  than  Tg.   Chapman  et  al  state 

C17,  pg  133: 

"the  apparent  radiant  temperature  is  much  closer 
to  the   flame   temperature  than  the  average   bulk 
gas  temperature. " 

Peak  flame  temperature,   Tr,  reached  in  high  r   CI  engines  is 

typically  4200  ±  200 ° R . 

3.2    KHAN  ET  AL  COMBUSTION  SPOT  PRODUCTION  MODEL 

While  a  valid  quantitative  calculation  of  soot 
concentration  in  New  Cycle  Engine  combustion  is  not  possible, 
a  brief  study  of  a  soot  production  model  for  a  simple  design, 
direct  injection,  quiescent  CI  engine  will  illustrate  the 
critical  parameters  on  which  soot  formation  is  based  and 
establish  a  foundation  for  empirical  estimation  of  peak  New 
Cycle    Engine    QR.     Khan   et   al   L181  have    developed 

semi -empiri cal  equations   for   soot  formation  in  a   fuel   jet 
mixing  model  applied  successfully  to  a  9.8  cm   bore  engine  at 
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FIGURE  4-1:  QR/QT  for  two  4.  stroke,  direct  injection  CI 

engines  under  varying  operating  conditions  as 
observed  by  Oguri  et  al  C83 
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2000   RPM.    This  simple  model  ignores  complex  stratification 

and   consideration  of  fuel  droplet  evaporation;   all   unburnt 

fuel  is  assumed  to  be  in  the  vapor  state.  The  jet  mixing  model 

injecting   into   a  cross-stream  air  flow  is  shown   in   Figure 

4-2.   The  air  entrainment  zone  (e)  is  composed  of  a  fuel  rich 

zone   (c)  where  soot  formation  occurs,   a  products   zone   (d) 

where   intense   reaction  occurs,   and  an  outer  zone   composed 

primarily   of  air.    Outside  the  entrainment  zone  is  the   air 

mass    zone    (a).     For   the   purpose   of   soot    formation 

calculations,   a   local  equivalence  ratio  of  unburnt  fuel   in 

the  fuel  rich  zone,   $  ,   is  defined.    A  series  of  equations 

are   proposed   by   Khan   »t   al   for   #u    calculation;    the 

significant    points   of   which   is   that   $u   is    inversely 

proportional    to    air   entrainment   which   is   in   turn    a 

progressive  function  of  air  swirl  ratio  and  engine  speed.   § 

is  not,   however,   directly  coupled  to  f  charge.   The  overall 

mechanism  of  soot  formation  under  high  temperature  and   short 

reaction   time   conditions  characteristic  of  diffusion   flame 

combustion    may   be   characterized   by   an   Arrhenius    type 

equation: 

d.S       V,,  -E=/RnT,  (52) 

u     -r        sou 

=  c  $  '-  P   e 

s      *u    u 

dt       V 

where:  dji  =  soot  formation  rate 
dt 

V  =  volume  of  soot  formation  zone  (c) 

V  =  cylinder  volume 

P   =  partial  pressure  of  unburnt  fuel  in  soot 

formation  zone  (c) 
T   =  local  temperature  of  soot  formation  zone  (c) 
c_  =  soot  formation  rate  coefficient 
E   =  soot  combustion  activation  energy 
R   =  universal  gas  constant 
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FIGURE  4-2:  Soot  production  fuel  jet  mixing  model  of 
Khan  et  al  D 8D 
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After  -formation,  soot  particles  coagulate  and/or  combust  from 
a  typical  initial  diameter  of  250  A  where  combustion  of  soot 
coagulate  is  again  controlled  by  an  Arrhenius  type  equation 
proportional  to  the  partial  pressure  of  0^. 

The   values   of   constants   and   coefficients  s.re,       of 
course,   a   function   of   individual  engine  design   and   fuel 
structure   selected   to  best  "curve  fit"   experimental   data. 
Valuable   qualitative  conclusions,   however,   stand  out   from 
the  work  by  Khan  et  al : 

1.  The  rate  of  soot  formation  depends  on  three 
principle  parameters:  local  P  (partial), 
T,  and  $  in  the  fuel  rich  soot  formation 
zone.  The  effect  of  varying  overall  § 
charge  is  moderated  as  §  is  always  near 
unity.  Swirl  ratio,  however,  has  a 
significant  effect.  A  comparison  of  soot 
concentrations  measured  by  Khan  et  al , 
Flynn  et  al ,  and  Dent  et  al  is  shown  in 
Figure  4-3  which  graphically  illustrates 
soot  concentration,  thus  luminous  flame 
radiation,  suppression  in  high  swirl  CI 
engines. 

Approximately  60"/.,  as  measured  by  Khan  et  al 
C18D,   of  soot  initially  formed  is  burnt  as 
it  moves  from  the  cooler  fuel  rich  zone  (c) 
into   the   products   zone   (d)   where   soot 
oxidation  incandescent   radiant  emission  is 
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a  grey  body,  optically  thick.  The  exhausted 
soot  is  that  portion  which  escapes  the  high 
temperature   products  reaction  zone  (d)   in 
the  process  of  turbulent  diffusion  from  the 
fuel  rich  to  the  lean  zone. 
3.     High   speed  photography  indicates  that  the 
majority   of  soot  formation  occurs   at   the 
beginning   of   the  main   combustion   phase, 
within    an   extremely   short   period   Csee 
Figure   4-33.    The   total  amount   of   soot 
remains   constant  in  the   later   combustion 
stages. 
A   clear   conclusion   of  this  model  is  that   both   soot 
formation   rate   and   soot  combustion   rate   are   progressive 
functions   of  fuel  and  0^  partial  pressures.    The  New   Cycle 
Engine,  operating  at  significantly  reduced  expansion  cylinder 
pressures,   should  therefore  experience  a  marked  reduction  in 
luminous  flame  radiant  heat  transfer.   Due  to  lack  of  (built) 
engine   test  data  however,   the  degree  to  which  QR   decreases 
relative  to  magnitudes  observed  in  high  r_  CI  engines  must  be 
estimated  from  indirect  empirical  relationships.   The  quality 
of   this   estimation,   based  principly   on   reduced   pressure 
consideration  in  low  swirl  conditions,   must  be  tempered  with 
the   knowledge   that  many  other  operating  parameters   have   a 
secondary  effect  on  flame  emissivity  and  temperature. 
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FIGURE  4.-3:  Comparison  of  in-cylinder  soot  conentration , 
£),  versus  crank  angle  as  observed  by  various 
investigators  0 OD 
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3 . 3    SECONDARY  FACTORS  AFFECTING  COMBUSTION  SOOT 

CONCENTRATION 

Extensive  tests  performed  by  Flynn  et  al  on  a  -Four 
stroke,  direct  injection  CI  engine,  with  piston  cup 
combustion  chamber  enlarged  to  suppress  swirl  ratio  and 
enhance  soot  production,  illuminate  many  of  these  secondary 
effects;  the  important  ones  of  which  a.re    listed  below  C163: 

1.  Qp  increased  with  inlet  manifold  pressure 
increase  at  nearly  proportional  rate; 
much  less  ignition  delay  and  fuel /air 
pre-mixing  occurs  causing  increased 
tendency  for  soot  formation  in  diffusion 
flame. 

2.  As  injection  time  was  advanced  an 
increase  in  Tg  and  an  increase  in  peak 
apparent  radiant  temperature  occurred; 
due  to  fourth  power  effect  a  significant 
increase  in  peak  radiant  emissive  power 
resul ts. 

3.  Decreasing  cetane  number  with  fixed 
reference  fuel  families  encouraged 
ignition  delay,  increased  fuel/air  pre- 
mixing,  and  reduced  soot  formation;  at 
supercharged  operation  this  effect  is 
effectively  suppressed  due  to  predominant 
effect  of  manifold  pressure  on  ignition 
del  ay. 
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4.    Increasing   C/H  mass  ratio  fuels  at  fixed 
cetan©  numbers  encouraged   Qp,  due   to 
simple      availability      of      carbon 
concentration;     the     rate    is     not 
proportional    to   C/H   mass   ratio    and 
apparent   only   at   high   ratio   aromatic 
fuel s. 
Flynn  et  al  additionally  note  that  the  apparent  flame  optical 
thickness,    krlp,    diminishes   approxi matel y    in    inverse 
proportion  to  increasing  cylinder  volume.   This  conforms  with 
the   Khan   et  al  model  and  Figure  4-3,   indicating   that   the 
number    of    incandescent   radiating   particles    in    later 
combustion  stages  is  essentially  constant. 
3.4    CHARACTERISTIC  HEAT  TRANSFER  EQUATIONS 

Emission   as   a   function  of  wavelength  for   clouds   of 
small  particles  obeys  a  variation  described  by  equation  (53): 

EA  =  1  -  exp(-k^lR/X'95)  (53) 

The  spectrum  is  near  grey,   and  for  simplicity  of  calculation 
can  be  converted  to  equivalent  grey  body  emissivity: 

r  n 
J*o 

e+  =  (54) 

qvSX 
Jxo 

where:     X   -  X   is  range  of  significant  X  band 

The  effect  of  combustion  cylinder  wall  reflectivity,   Tw?   as 

demonstrated   by  Oguri  et  al  in  Figure  4-1,   may  be  corrected 
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-for   by  determination  of  an  apparent   -flame   emissivity,   £_, 
according  to  equation  (55): 


ci-Ve*       "Mr 

Ea  = =   1  -  e  (55) 

^w^f 

where:      k^   =  apparent  flame  absorption   coefficient 

per  unit  length 

When  -flame  envelopes  the  entire  cylinder  volume,   observed  by 

Alcock   et   al   CI 13  to  occur  approximately   10  to  20°   after 

ignition   start,   1 R  may  be  calculated  by  equation  (47)  as   a 

function  of   instantaneous  cylinder  V  and   Amt.   Sitkei  et  al 

C15H  state  that  surfaces  of  pre-swirl  chambers,   piston  crown 

and   cylinder   heads  have  ?w  =  0.1  (appro::.  );  cylinder  sleeve 

^   =  0.2  (appro:-!.  )  .   It  should  be  noted  that  estimates  for  T 

vary   greatly,   in  part  due  to  different  modeling  of  cylinder 

surfaces   as   specular   reflective   to   semi -rough   diffusive 

reflective   characteristic  of  oil   soaked,   soot-contaminated 

surfaces.   Qr,  may  then  be  determined  by  equation  (56): 

QR  =    eacrQAHTCT+4-Tw43  <1    hr/3600    sec)  (56) 

^BTU/sec 

3.5    SITKEI  ET  AL  METHOD  OF  APPARENT  FLAME  EMISSIVITY,  E_, 

CALCULATION 

A  useful  empirical  relationship  used  by  Annand  [3]  and 
Dent  et  al  C103  based  on  known  soot  concentration  is  stated 
in  equation  (57)  in  U.S.  units: 

kf  =  5898  &    ft-1  (57) 

where:  &    =  soot  concentrations  in  lbm/ft'"' 
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For  estimation  of  £a  without  knowledge  of  soot 
concentration,  Sitkei  et  al  propose  the  following  definition 
of  £a  as  a  function  of  instantaneous  P,  $,  1 ,-.  and  crank  angle 
past  ignition  start  C153: 


-.  3048  kF'lr-, 
P   R. 


E,  =  £_  CI  -  e  3  (58) 

a     o 

where:  £Q  =  emissivity  of  infinitely  thick  radiating 
f  1  ame ;  ss  0 .  8 

k    =   apparent  flame  absorption  coefficient   per   unit 
length-unit  pressure  (atm) 

Values  of  k  Are  determined  from  Figure  4-4,  at  applicable 
crank  angle,  and  excess  air  factor  equivalent  to  §  .  Sitkei 
et  al  based  Figure  4-4  on  measurements  conducted  by  a  photo- 
electric pyrometer  unit  through  a  quartz  view  window  insert 
to  a  swirl  pre-combustion  chamber  of  a  four  stroke  Steyr 
diesel  at  r  equal  to  18.5.  Recognizing  that  £_  is  dependant 
on  many  additional  variations  in  engine  operating  parameters, 
Sitkei  states  that  to  a  first  approximation  C14,  pg.  563: 


"in  the  absence  of  other  data,  it  is  felt  that  k 
given   here   is   universally   adaptable   to    all 
engines  of  the  diesel  type....  with  the  known  P-V 
diagram  of  the  engine  under  test  the  variation  of 
k   can  be  found" 

Inherent   in  this  assumption  is  an  assumed  %       of  10"/.  to   207., 

an  ignition  start  near  TDC,   a  rather  arbitrary  value  of  £   = 

0.8   based   on   observed   steady   state    diesel   fuel   flame 

emissivities   in  combustors,   and  the  belief  that  £   observed 

in   a   2.6  cm  pre-combustion  chamber  may  be   extrapolated   to 

main  cylinder  emissivities. 
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+  0  from  ignition  start  I     decrees 


FIGURE  4-4-1  Variation  of  kD  versus  angular  displacement  for 
Ca  calculation";  Sitkei  et  al  model  L"KJ 
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It  should  be  noted  that  the  values  of  E   determined   by 

Cm 

the  above  method  are  generally  less  than  observed  by  a  number 
of  other  investigators  in  direct  injection  CI  engine  testing. 
Dent  et  al  state  II 10,  p.  1760  3: 

"The  essential  results  of  Sitkei  and  Ramanaiah 
are  in  agreement  with  those  of  Flynn  et  al  in 
terms  of  the  trend  of  flame  emissivity  and 
optical  thickness  with  crank  angle,  load,  and 
overall  air  fuel  ratios.  However,  the  magnitude 
of  the  flame  emissivity  and  optical  thickness  ar& 
approximately  half." 

For  first  approximation  of  predicted  peak  E_  in  the  New  Cycle 

Engine  expansion  cylinder,   it  is  therefore  believed  that  the 

E_   obtained  by  the  Sitkei  et  al  method  should  be   multiplied 

by   a  factor  of  2.    It  is  also  noted  that  Sitkei  et  al  agree 

with   typical   QR/Qy  range  of  207.  to  30"/.  even  at   reduced   £ 

i 
values  due  to  assumption  of  Belinskiy's  significantly   higher 

Tr    profile   shown   in   Figure   4-5,    vice   direct   optical 

measurement . 

3 . 6    LQUJ  COMBUSTION  PRESSURE  POWER  UNIT  METHOD  OF 

APPARENT  FLAME  EMISSIVITY.  £  ,  CALCULATION 

I—  .1.—  I.. -.■—.■■■■  —   ■ — M ■ i  — ...  ■   in.  ■      I ^   f.   ■   II.. — — ■■■■■...l 

For  first  approximation  purposes,  an  alternate  method 
of  predicting  £,  is  to  examine  the  characteristic  emissivity 
equations  used  in  power  units  operating  at  the  low  peak 
combustion  pressures  of  the  New  Cycle  Engine.  Different 
properties,  continuous  flaming  and  non-reciprocating 
mechanics  of  gas  turbines  do  not  allow  broad  application  of 
flame  emission  data  to  internal  combustion  engines.  However, 
actual  combustion  kinematics  and  flame  luminosity  observed  in 
power  units  at  pressure  conditions  more  reflective  of  the  New 
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FIGURE   1-5: 


Measured  Tf  in  CI  engines  versus  crank  angle  as 
observed  bv  various  investigators  p 0 J 
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Cycle   Engine   should   not   be   ignored.     It   is   therefore 
considered  that  a  fair  prediction  of  New  Cycle  Engine  peak  £ 
only   can   be   ascertained  by  use  of  gas  turbine   power   unit 
emissivity  equations  using  diesel  fuel  properties. 

For  gas  turbine  application,  £  can  be  determined  by 
equation  (59)  CI 9 3: 

£=  =  1  -  e~z  (59) 

where:   2  =  -33636.55  PC  (fuel /air  mass  ratio)  1R3  '  5Tg~1 '  5£ 

The  flame  luminosity,  £,  can  be  determined  as  a  function  of 
fuel  C/H  mass  ratio  from  Figure  4-6.  For  commercial  liquid 
fuels  of  conventional  character  ranging  from  gasoline  to  fuel 
oil,  C/H  mass  ratio  can  be  estimated  using  equation  (60): 

C  mass  ratio  =  .  16015/>  -  1.66  (60) 

H 

where:   p  fuel  density  in  lbm/ff"'  is  referenced  to 
518° R  (288° K) 

3.7    BOLTZMANN  NUMBER  FUNCTION  METHOD  QF  FLAME 

Temperature,  Tf ,  CALCULATION 

A  consensus  opinion  of  investigators  support  that  for 
first  approximation  purposes,  radiant  temperature  is  equal  to 
Trj  various  Tr  profiles  published  from  high  rc  CI  engine 
testing  are  shown  in  Figure  4-5.  An  empirical  correlation 
allowing  calculation  of  T^  is  the  Boltzmann  number  function, 
Tg/Tf  =  f  (Bo#) ,  established  in  Figure  4-7.  The  Bo#  is 
defined  in  accordance  with  equation  (61)  C143: 


86 


^ 


2  4  6 

fuel    C/H      TTidss  ratio 


y 
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±t      BeCp  fuel 
Bo*  =  _  (61) 

■-xoef  AHTTf 
where:     Be  =  fuel  consumption  rate 

C   fLiei  =  mean  total  specific  heat  of  combustion 
products  of  1  lbm  of  fuel 

For  ease  of  calculation,  Be  and  C   fuel  may  be  calculated  for 

the  New  Cycle  Engine  using  the  following  identities: 


(SFC) (Aheat)  <mcharqe>  (Nrpm>  (60  min/hr> 

Be  =  — ° (62) 

lbm/hr  fuel  2544  BTU/hr-HP 


C„  fuel  a;  (Air/Fuel  mass  ratio)  (C„  combustion  products)   (63) 

Substitution   of   the   above  identities   into   equation   (61) 

yiel ds: 

(64) 

C(SFC)  (Aheat)  (m.=  „__)  <NC.C.M)  (Air  /fuel    mass    ratio)C 
+t                                                   charqe         Kr-rl  p,-n 

Bow    =   : up 


£f AHTTf " 


where:  C  =  1.3768xl07 


—  # 

Calculations  of  T,  at  a  known  lg  and  £r  from  the  Bo   function 

requires   an   iterative  process  matching  an  assumed   T,   both 


directly  from  the  generated  Tg/Tr  fraction  and  its  coupled 
Bo  from  Figure  4-7.  An  increasing  £r  will  decrease  the  Bo  , 
thus  decrease  Tg/Tr,  resulting  in  an  increased  calculated  Tr. 

It   is  noted,   however,   that  increases  in  E^  and  T^  are   not 

+t  _ 

proportional;  over  a  wide  range  of  Bo   (  >  .  6  ),  Tg/T^  slope 

is   only   slightly   positive.     At   Bo    greater   than   1.2, 

characteristic   of   low  luminous  flames,   the  value  of  T^   is 

i±  — 

essentially  independent  of  Bo   (thus  £+)  and  equal  to  Tg/0.8. 


S3 


Bo      number 


FIGURE  4-7:  Variation  of  Tg/Tf  in  the  function  of  the  Bo 
number  {J4-J 
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3.S    PREDICTION  QF  PEAK  COMBUSTION  S=  and  T, 

— — ■■  ■  -——--■■-— — — — — — ■■— — —c\  '■■■ -■■  —■■■■■*!* 

For   calculation   o-f   peak  Ea  and  Tr  in  the  New   Cycle 

Engine   expansion   cylinder   at   full   load,    the  following 
assumptions  are  made: 

1.  The  expansion  cylinder  reflectivity,  %  =  0.1 

vv 

2.  The  flame  growth  encompasses  the  entire 
cylinder  volume  at  20°  maximum  crank  angle  from 
ignition  start  Ccrank  angle  =  65° J;  flame 
growth  rate  is  constant. 

3.  The  flame  luminosity,  £,  calculation  is  based 
on  medium  speed  diesel  fuel  C/H  mass  ratio  =  6.3 

4.  The  fuel /air  mass  ratio  =  .05  (1/20)  for  $  =  . 75 

5.  C„     at   fuel   C/H  mass  ratio   =   6.3   mav   be 

Pep 

approximated   by  equation  (32)  as  a  function  of  Tg 

for  CnH^n  fuels  (C/H  mass  ratio  =  6.0). 

The  Sitkei  et  al  method  of  £_  calculation  was  developed  on   a 

test   engine   exhibiting  typical  high  r   combustion   pressure 

profiles,  and  TDC  ignition  start,  both   conditions  not  met  in 

New  Cycle  Engine  operation.    The  peak  k   from  Figure  4-4   at 

the  correct  excess  air  factor  value  will  therefore  be  applied 

at  the  crank  angle  of  maximum  combustion  pressure  where  flame 

growth  reaches  1007.  cylinder  volume. 

Sitkei  Method;  £_  peak 

8r      ,  =  65° 
£a  peak    =^=— 


excess  air  factor  =  (.75)    =  1 . ZZ 

peak  k   at  1.33  excess  air  factor  =  0. 63 

(see  Figure  4-4) 
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Using    equations     (47)     and     (58),       plus    applying    a    £       correction 
•factor    of    2: 


8  P  V  AHT  1R  kp  .3048kpPlR       £a  2E, 

ati  ft3         ft2  ft     (ft-at*)"1 

65'    3.76064      .488331      3.944       .4210         .63  .304  .210       .420 

peak    £a    =    .420 
at    8    =    65° 


Low    Combustion    Pressure    Power    Unit    Method,    £       peak 

■flame    luminosity,    <£    =    4. 0 
(see    Figure    4-6) 

Using    Equations     (47)     and     (59): 

9           P  Tg               V             A       flaie  growth      lRiX  Z          Ea 

ati  »R             ft3  it2  I  H 

SO  4.24684  22!?. 6  .309053  3.358  25  .0782  .382  .318 

55  4.09172  2571.21  .366348  3.545  50  .1757  .455  .366 

60  3.92401  2864.13  .426316  3.741  75  .2906  .477  .379 

65  3.76064  3136.02  .488331  3.944  100  .4210  .480  .381 

70  3.26107  3064.15  .551767  4.152  100  .4518  .447  .360 

peak    £a    =    .381 
at    6    =    65° 

A   repeated   calculation  -for  both  methods  at  peak   combustion 

pressure  greater  than  50  atm,   typical  -for  direct   injection, 

high  rc  CI  engines,   results  in  peak  £a  m    1.0,   which  closely 

agrees  with  the  expected  value  (:s:0.9). 

Boltzmann  Number  Function  Method,  peak  T^ 

The   known  values  required  -for  calculation  o-f  the  Bo   at   65° 

crank  angle,  peak  £_  o-f  (.420  +  .38D/2  =  .40  Are: 
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Tg  ■  3 1 36 .0  °  R 

SFC  =  .284298  lbm  -fuel/HP-hr 

Aheat  =  500  BTU/lbm  charge 


4-    mcharge  =  ■ °24537  lbm 


NRPM  ~  813-3 

C„     =  . 247+.243x 10~4 (3136.0° R)  =  .3232  BTU 


Pep 


lbm-°R 

7.  air  /-fuel  mass  ratio  =  20 

8.  AHJ  =  3.944  -ft2 

Substituting  the  above  values  into  equation  (64)  yields: 


Bo#  = 


(1.3768xl07) (.28429) (500) (.024537) (818.3) (20) (.3232) 

(3.944)£.fTf3 


6.435x 1010 

«  Y    3 

b-f  '-f 

From   equation   (55)  ,   an   estimated  peak  £_  o-f   0.40   at   an 

CI 

assumed   ?    =   0.1  yields  a  peak  C^   equal   to   0.426.    The 
predicted   peak  Tr  may  there-fore  be  calculated  -from  iteration 

tk 

using  Figure  4—7  at  Bo   equal  to: 


Bo#  = 


6.435xl010  =  l.SllxlO11 


.426Tf-         Tf° 


Assuming  a  value  of  T+  =  Tg/0.8  =  3920° R,  the  calculated  Bo# 
equal  to  2.51.  This  Bo  is  greater  than  1.2,  there-fore  the 
estimated  peak  T^  equals  3920°R.  Rounding  o-f-f  -for  -first 
approximation: 

peak  T,  =  4000° R 
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To 


check   the  accuracy  of  this  procedure,   a   peak   £ r 


value  of  1.0  is  assumed  commensurate  with  observed  peak  E_  in 

d 

direct  injection,  high  r   CI  engines: 


Bo*  = 


010  =  6.43.-" 


6.435x 10* "  =  6. 435x 10 


(1.0)Tf-        Tf° 


From  Figure  4-5,   the  observed  peak  T,  is  typically  equal   to 
4200° R;   Tg/Tf  now  equals  0.75.   Substituting  4200° R  into  B~ 


o 


=  6.435X  1010/Tr3  yields  a  value  o-f  0.87  which,  entering 
Figure  4-7,  extrapolates  to  Tg/Tr  =  0.76.  The  iteration 
essentially  matches,  indicating  that  a  peak  £_  and  peak  Tr 
equal  to  1.0  and  4200°R  respectively,  are  in  fact  reasonable 
values  -for  direct  injection,  high  rc  CI  engines. 
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4.      CALCULATION  QF  ORDER  QF  MAGNITUDE  REDUCTION 

IN  COMBUSTION  RADIANT  HEAT  FLUX 

Woschni  did  not  separately  measure  radiant  heat  flux  in 
his  engine  testing;  for  ease  of  calculation  he  simply 
absorbed  the  radiant  component  into  the  forced  convection 
heat  transfer  model  by  use  of  time  averaged,  "curve  fitting" 
multiplicity  coefficients.  Woschni  in  fact  believed  that 
C9,  pp.  3074  -  30753: 

"the  total  radiation  amounts  to  no  more  than  207. 
of  the  total  combustion  term. .. In  steady  flames 
there  are  usually  large  combustion  chambers  with 
great  volumes  of  flames,  and  because  of  the 
small  velocities  there  are  small  convective  heat 
transfer  rates  only,  so  that  in  these  cases 
radiation  is  responsible  for  an  appreciable 
amount  of  the  total  heat.  Because  of  the  high 
velocity  and  density  of  gas,  the  convective  heat 
transfer  coefficients  in  internal  combustion 
engines  are  about  ten  times  higher,  so  that 
despite  the  higher  temperatures  the  radiation 
plays  an  appreciably  smaller  part." 

While  the  thrust  of  the  statement  is  correct,  Woschni 
unfortunately  underestimated  Qn,  basing  cursory  radiant  heat 
transfer  calculations  on  the  assumption  of  £_  =  0.6  from 
steady  state  flame  observations  and  employing  for  the  unknown 
radiant  temperature,  Tg ,  vice  T^ . 

For  the  purpose  of  this  calculation,  establishing 
Woshni  peak  £a  and  T^  values  of  1.0  and  4200° R  better 
conforms  to  published  data  on  CI  engine  testing  of  similar 
design,  direct  injection,  high  r  characteristics.  It  is 
noted  that  actual  observed  peak  £_  values  are  typically  equal 
to  0.9  vice  1.0.  To  maintain  continuity  with  £_  prediction 
methods   used  in  section  3.8,   however,   a  value  of  unity   is 
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assumed   as   predicted   -For  high   peak   pressure   conditions. 
Without   knowledge  o-f  the  radiant  component  o-f  qw  and  q^ipp  as 
a   -function  o-f  crank  angle,   it  is  not  practical   to   attempt 
calculation  o-f  a  mean  order   o-f   magnitude   reduction  ratio. 


Cu^^NCE '  over  the  effective  90°  crank  angle  domain  where 
luminous  radiant  heat  transfer  can  be  significant  (£^  ^  0.1) 
in  high  r  ,  medium  speed  CI  engines.  Fortunately,  in  this 
case  it  is  not  required.  Making  the  reasonable  assumption  of 
equivalent  £_,  Tr  decay  rates,  flame  emissivity  in  the  New 
Cycle  Engine  expansion  cylinder  should  certainly  decay  below 
0.1  at  advancing  crank  angles  from  ignition  start  prior  to 
+90°.  Thus,  peak  qij/qMrF  should  be  a  mini  mum ,  conservative 
estimation  of  the  mean.   Therefore,  from  equation  (56): 

ea<T/-Tw4)w 


qw 

« 

qW 

qNCE 

qNCE 

rad 

L  at 

ion 

peak   « 

ea<T/-Tw4)NCE 


T,  =  900° R 
w 


peak 


Substituting   in   the  respective  Z&    and  T^  values   previously 
determined: 

qw         1.0(42004-9004) 

minimum  = —   =   3.047  approx . 

qNCE      0 . 4 ( 4000-900  ) 


CHAPTER  FIVE 


CALCULATION  OF  CORRECTED  HEAT  TRANSFER  LQSS/HEAT 

RELEASE  PER  CYCLE,  ft* 


Total  heat  transfer  loss,  Qj,  tor  a  natural 
aspirating,  four  stroke  CI  engine  includes  heat  transfer 
during  the  induction  and  compression  strokes  which  do  not 
exist  in  the  two  stroke  New  Cycle  Engine  expansion  cylinder 
cycle.  In  the  induction  stroke,  the  heat  transfer  driving 
potential,  Tg-Tw,  is  negative;  induction  charge  heating  is 
boosted  by  residual  gas  mixing.  In  the  compression  stroke, 
adiabatic  compression  may  boost  Tg  to  greater  than  1500° R. 
From  Figure  2-1,  it  is  observed  that  a  significant  positive 
Tg-Tw  occurs  only  in  the  last  30°  of  the  compression 
stroke.  The  order  of  magnitude  of  this  positive  driving 
potential  compared  to  that  observed  in  the  combustion- 
expansion  and  exhaust  strokes  is  still  guite  small.  It  may 
therefore  be  assumed  that  the  induction  and  compression 
stroke  contribution  to  Qj  for  a  natural  aspirating,  four 
stroke,  high  r  CI  engine  is  not  significant.  A  direct 
coupling  of  observed  QR/Q-j-  ratios  may  then  be  applied  to  the 
New  Cycle  Engine  expansion  cylinder. 

From   Chapters  2,   3,   and  4,   the  following   principal 
values  have  been  determined: 


1.  ftT  uncorrected  =  19.7617.;  Woschni  correlations 

2.  QR/QT  Woschni  engine  ^  V;  .  2  ^  V  ±=  .  3 
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Forced  convection  q^/c,NCE  =  6,<-)87 


4.  Radiation    qui/dure-    =    3.047    mi  ni  mum 

fly  must  be  corrected  tor  the  -failure  of  Woschni  's 
correlation,  based  solely  on  a  theoretical  forced  convective 
heat         transfer       model,  to      properly       scale       down       luminous 

radiation    heat    transfer    to    low    r       operating    conditions. 

Correction    to    Tlj    =    Ttpp    +    Tip 
At    high    r       operation: 


Qc,  Qc./  (Aheat )  (m    .   ____)  Tlr-. 

K  K  charge  K 


=    V 


QT  <QR    +    QFC)/(^heat)(mcharge)  TtR    +    1lpc 

HR    =    WTlR    +    WRFC 

TlR(l-¥)     =    VTlFC  ' 

HR    =    _W_   TlFC  (65) 

1-V 

Substituting  into  equation  (65)  the  known  order  of   magnitude 

reduction    ratios   for   low   r    operating   conditions    and 

establishing  equality: 

At  low  r   operation: 

3.047  IU*  =  _¥_  6.087  tlFC* 
1-W 


TlR*       =  2.00   V   ftFC*  (66) 

maximum         1-W 


As    Woschni 's    correlation   properly   scales   down    forced 
convection  heat  transfer  only,  ftFC*  =  (1  -  W)TlT.   Therefore  a 
conservative,   first   approximation  prediction   of   corrected 
total   heat   transfer   loss/heat   release   per   cycle,   Tlj  , 
for  .2  £  W  i  .3  is  : 
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Forced  convection  heat  transfer  loss,  UFC* 
HFC*  =  <l-¥)1lT  =  <1-W)  19.7617. 
13. 857-/.  £  UFC*  £     15.8097. 
Radiant  heat  transfer  loss,  TtR 

TtR*  =  2.00(  V  >%c*  =  2,00<  ¥  )  (l-V)tlT 

1-W  1-V 

=  2 . 00  ¥  HT 

=  2.00  W  19.7617. 

7.9047.  6  Up*  g  11.8577. 

max  i  mum 

Total  heat  transfer  loss,  tlj 

RT*  =  ftR*  +  tlFC*  =  2.00WR-J-  +  <l-¥)flT 

=  (2.00V+1-W)RT 

=  (l+W)TlT 

=  (l+¥)  19.7617. 
23.7137.  A    g.T*  g  25.6897. 


max l mum 
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CHAPTER  SIX 

CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  STUDY 

A  semi-empirical  correlation  -for  mean  instantaneous 
heat  transfer  coefficient,  developed  by  Woschni  on 
conventional  four  stroke,  direct  injection,  quiescent  CI 
engines,  was  modified  to  predict  heat  transfer  loss  in  a 
computer-generated  preliminary  design  engine  cycle.  The 
engine  divides  the  four  stroke  functions  between  a 
compression  cylinder  and  an  expansion  cylinder  interconnected 
by  a  static  regenerative  heat  exchanger.  Due  to  regenerator 
boost  heating  of  the  air  charge,  required  compression  ratios 
in  the  New  Cycle  Engine  expansion  cylinder  will  be 
significantly  reduced,  and  increased  thermal  efficiency 
achieved.  To  verify  applicability  of  Woschni 's  correlation 
to  expansion  cylinder  heat  transfer,  points  in  the  cycle  were 
examined  comparing  Woschni 's  correlation  to  fully  developed, 
turbulent  flow,  forced  convection  relationships  for  flat 
plate  and  pipe  flow.  The  mean  instantaneous  gas  velocity  was 
estimated  from  the  construction  of  the  charge  impulse 
transient  existing  into  the  expansion  cylinder  from  the 
regenerator  port.  An  order  of  magnitude  reduction  in  forced 
convection  and  luminous  radiant  heat  flux  was  calculated  to 
investigate  the  ability  of  Woschni 's  correlation  to  scale 
down  correctly  total  heat  transfer  to  low  compression  ratio 
operation. 
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It  was  found  that  the  impulse  velocity  transient  was 
suff  i  ci  entl  y  suppressed  to  approximate  quiescent  conditions. 
It  was  -found  that  laws  of  similarity  incorporated  into 
Woschni's  correlation  appear  to  allow  uncorrected  application 
at  low  compression  ratio  operation  -for  -forced  convection  heat 
transfer  only.  The  order  of  magnitude  reduction  in 
combustion-expansion  stroke  forced  convection  heat  flux  was 
found  to  be  as  much  as  twice  the  order  of  magnitude  reduction 
in  radiant  heat  flux.  A  correction  was  proposed  for  the 
application  of  Woschni's  correlation  to  achieve  a  conservative 
prediction  of  heat  transfer  loss  in  the  New  Cycle  engine 
expansion  cylinder  cycle: 


QT  =  1.25  QT  (67) 

corrected  uncorrected 

where:  QT  =  f  Q  d© 

uncorrected    J 

Q  =  he  AHT  (Tg-Tw) 

he  =  15.480  P*8  CC^Sp  IT8  B~'2  Tg~*53 

C*     =  6.18-      exhaust  stroke 

2.28      combustion-expansion  stroke 


This  study  investigated  heat  transfer  in  engine 
operating  conditions  significantly  below  the  range  for  which 
published  correlations  exist.  Theoretical  work  and  actual 
measurement  of  an  engine  prototype  remain  to  be  completed  to 
prove  the  validity  of  significant  extrapolation  of 
reciprocating  engine  heat  transfer  correlations.  This  study 
points  out  three  specific  areas  requiring  future  study 
concerning  reciprocating  engine  heat  transfer  prediction: 
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1.  Quantitative  prediction  of  -flow  pattern  magnitude 
and  distribution  as  a  -function  of  crank  angle  and 
1 ocation. 

2.  Quantitative   prediction  o-f  soot  concentration, 
apparent    flame    emissivity,     and     radiant 
temperature    as    a   function   of    operational 
variables    in   low   pressure   diffusion    flame 
combustion. 

3.  Quantitative    effect   of    appreciable    fuel 
injection  and  ignition  delay  in  the  expansion 
stroke  of  low  compression  ratio  reciprocating 
engines  on  incandescent  particle  concentration 
and  luminous  flame  radiation. 
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APPENDIX  A 

The  purpose  of  this  appendix  is  to  present  preliminary 
design  parameters  and  the  computer-generated  thermodynamic 
■firing  cycle  -for  the  New  Cycle  Engine  expansion  cylinder  on 
which  the  heat  transfer  calculations  are  based  Zll: 

Preliminary  Design  Parameters;  expansion  cylinder 

1.  Mean  piston  speed,  3   =  2000  ft/min  (33.33  -ft /sec) 

2.  Square  con-figuration;  stroke  =  bore  diameter  (L=B) 

3.  Neglect  clearance  volume;  V'c  =  0 

4.  Neglect  cylinder  sleeve  clearance 

5.  Regenerator   port  exit  cross-section  area,   A,-.  = 
0.80  ft   (l/3rd  cylinder  circumference  x  .75  inch 
depth) 

6.  Compression  cylinder  supercharge  pressure  ratio=l; 
natural  aspirating 

7.  Quiescent   -fuel   injection  at  nominal   25°   crank 
angl  e 

8.  Ignition  start  at  nominal  45°  crank  angle 

9.  Specific   heat  release  per  cycle,  Aheat       =   500 
BTU/lbm  charge. 

10.  Specific  fuel  consumption,  SFC  =  .284298  lbm/HP-hr 

11.  Fuel  base  =  nominal  C  H-,   medium  speed  diesel  fuel 

12.  Full  load  equivalence  ratio,  $  =  .75 

13.  Full  load  air/fuel  mass  ratio  =  20 

14.  Mean  heat  transfer  surface  wall  temperature, 
Tw  =  900° R 

15.  Specific    gas   constant   for   air   charge    and 
combustion  products,  R  =  53.90  1 bf -f t/lbm-° R 

16.  Air  charge  Prandtl  number,  Pr  =  .70 
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TABLE  A-l 


New   Cycle  Engine  expansion  cylinder  computer-generated 
thermodynamic   -firing   cycle   model  as  a   -function   o-f   crank 


angle;  5°  increments  CI] 


crank  angle 


Pressure 
P,  atm 


mean  gas  temp 


Vol ume 


X 


ft- 


EXPANSION-INDUCTION  CONDITIONS 


10 
15 
20 
25 
30 
35 
40 
45 


2.49823 
3. 1201 
3.91693 
4.83479 
5.67408 
6.10171 
5.88219 
5. 12078 
4. 14883 


2149.04 
2159.25 
2195.73 

nnne   *-*  / 
*   ■*'  J  ^7  •  ilO 

'-»'->'->0   "7  1 

jL.jL.jL-V  .  /  1 

2194.29 
2116.61 
2010.57 
1896.6 


3.4067E-0; 
. 0 1 35883 
. 0304243 
.0537145 
.0831858 
. 118493 
. 159226 
.204916 
. 255045 


EXPANSION-COMBUSTION  CONDITIONS 


50 

4.24684 

2251.6 

. 309053 

55 

4.09172 

2571.21 

. 366348 

60 

3.92401 

2864. 13 

.426316 

65 

3.76064 

3136.02 

.488331 

70 

3.26107 

3064. 15 

.551767 

75 

2.88976 

2998.73 

. 616005 

80 

jl.  ■  jd66i! 

2939. 11 

. 680446 

85 

2.33694 

2884.75 

.744518 

90 

2. 12966 

2835. 18 

.887682 

95 

1.95638 

2789.98 

.869441 

100 

1.81064 

2748.78 

. 929342 

105 

1.68741 

2711.25 

.98698 

110 

1.58277 

2677 . 09 

1 .  042 

115 

1.49361 

2646.03 

1 . 09408 

120 

1.41745 

2617.83 

1. 14298 

125 

1   TCOTO 
1  •  OJ<lOZ 

2592.29 

1. 18847 

130 

1.29661 

2569.22 

1 . 23038 

135 

1 . 24904 

2548.44 

1.26857 

140 

1 . 20857 

jL.*_l  JL.  7  .  QL 

1.30291 

145 

1. 17433 

2513.21 

1 . 33334 

150 

1. 14565 

2498.51 

1 . 3598 

155 

1. 12197 

2485.62 

1 . 38223 

160 

1. 10286 

2474.46 

1 . 4006 1 

165 

1 . 08796 

2464.96 

1.41492 

170 

1 . 07702 

2457.07 

1.42515 

175 

1 . 06986 

2450.74 

1.43129 

BDC 

1 . 06636 

2445.96 

1 . 43333 
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TABLE  A-l  (CONT'D) 


8°        P,  atm  Tg,  °R  V,  ff' 

EXHAUST  CONDITIONS 


BDC 

1.05 

185 

1 .  05 

190 

1 .  05 

195 

1.05 

200 

1 .  05 

205 

1.05 

210 

1 .  05 

215 

1.05 

220 

1 .  05 

1 .  05 

230 

1 .  05 

235 

1 .  05 

240 

1 .  05 

245 

1.05 

oc;a 

1.05 

255 

1 .  05 

260 

1 .  05 

265 

1.05 

270 

1 .  05 

275 

1 .  05 

230 

1.05 

285 

1 .  05 

290 

1 .  05 

300 

1 .  05 

305 

1 .  05 

310 

1 .  05 

315 

1.05 

320 

1 .  05 

325 

1 .  05 

330 

1 .  05 

335 

1 .  05 

340 

1 .  05 

345 

1.05 

350 

1 .  05 

355 

1 .  05 

TDC 

1 .  05 

2439.06 

1 . 43333 

2434.63 

1.43129 

2430.7 

1.42515 

2426.81 

1.41492 

2422.95 

1.40061 

2419. 13 

1 . 38223 

2415.36 

1 . 3598 

2411.64 

1 . 33335 

2407.99 

1 . 30292 

2404.41 

1.26858 

2400.9 

1 . 23039 

2397.47 

1. 18849 

2394. 14 

1.  143 

2390.89 

1.0941 

2387.75 

1.04201 

2384.71 

.985995 

2381.77 

. 929358 

2378.95 

. 863458 

2376.24 

.  8077 

2373.65 

. 744536 

2371. 18 

. 680464 

2368.82 

.616023 

2366.58 

.551785 

2364.46 

. 488348 

2362.44 

. 426333 

2360.54 

. 366364 

2358.73 

. 309069 

2357 . 03 

. 255059 

2355.41 

. 204929 

2353.88 

. 159238 

2352.41 

. 118504 

2351.01 

.083195 

2349.67 

. 053722 

2348.36 

. 0304295 

2347.09 

.0135921 

2345.83 

3.4083E-0 

2344.58 

0 

io: 


0  I   dear 


grees 


FIGURE  A-1 :  Generation  of  expansion  cylinder  firing  cycle 


motoring  P0 
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APPENDIX  B 

The    purpose   o-f   this   Appendix   is   to   present   the 
derivation   of   several   geometric   relations   -for   a   simple 
geometry   reciprocating  cylinder.    Figure  B-l  illustrates   a 
standard  cylinder  with  de-fined  piston  linkage  symbology. 
1 .  Instantaneous  piston  speed,  Sp 
by  the  law  o-f  cosines: 

lr2  =  s2  +  a^  +  2ascos6 

s^    -  (2acos8)s  +  (a2-l  2)  =  0 
by  the  quadratic  -formula,  solving  -for  s: 

-<-2acose)  ±  C(2acose)2  -  4(a2-l  2)21/2 


2 

s  =  acose  ±  Clr2  +  a2(cos2e-i) ]1/2 

substituting   in  the  trigonometric   identity,   cos"^6-l=sin^e, 
and  taking  the  positive  real  root: 

s  =  acose  +  Clr2  -  a2sin26Il1/2  (B-l) 

di -f -f  erenti  ati  ng  s  with  respect  to  time,  t: 

ds  d8     #  ds 

dO  dt       d8 


2  •  -    - 

a  sinycosy 

Sp  =  eC-asinO  -  3 

<lr^-a^sin^8) l/^ 

substituting  in  the  crank  ratio  identity,  Rc  =  lr/a,  and 
simpl i  f ying: 

cose 
Sp  =  -aesineci  + 


"7  "7        1  /' 

(R^-sin^e)  x/ ■ 


substituting    in    the    e    identity,    e    =    2irNRpri/60: 
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-aZirNpp^sinU  cosO 


Sp    =       C     1    +    1  (B-2) 

60 

2.     Mean    Piston    Speed,     Sp 


60  (R    *-ain2e)  1/2 


rir 

Sp   s    l_l    Sp    de 

IT     J 


o 


1       ^      -a2irNRpMsinO  cosO 

=  _  F  C1+ 3  dO 

0        60         (Rc^-sin^0)  i/':- 


-2aNr-.c,M  _  _     sinOcosO 

Rrn  pir  rir 

C  sinO  dO  +  - - T^e    1 

60  JQ  J0  (Rc2-sin20) 1/2 


■2aNr,r.M         «-     sinOcosO 

[  2  +     d8  3 


7      "?         1  /*? 
60  JQ  (Rc^-sin^O) l/^ 


Let  u  =  cosO  and  substitute  into  the  integrand  the  equivalent 
identx ti  es: 

du  =  -sine 

sin^e  =  1-cosO  =  1-u"1" 

cos(tt)  =  -1 

cos(O)  =  1 


-2aNRpM  -        U 

Sp  =  C  2  +  I    du  : 

60  J+l  (R  ^-l+u*-)  i/j:- 


-2aNRpM  t 

C  2  +  I    du  3 

60         J_1  (Rc^-l+ur) l/- 


(  2  +  CRc^--i+u^:i/^  ) 

60  - 1 
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2alMRpM 

Sp  =  (2+0) 

60 


Sp  = 


-4aNRpM 


60 
Substituting  in  the  crank  radius  identity,  a  =  L/2: 

2LNRPM 
Sp  =  (B-3) 

60 

3 "  Instantaneous  Cylinder  Volume,   V : 

Neglecting  clearance  volume,  V  ,  and  sleeve  clearance: 

V  =   ttB2  Le  (B-4) 

4 

4*  Instantaneous  Heat  Transfer  Area,   A|_iT: 

Neglecting  clearance  volume,  Vc ,  and  sleeve  clearance: 

AHT  =  2<H*2)  +  nBLe  (B-5) 

4 

where:  L^  =  1 _  +  a  -  s 
t»     r 

s  -from  equation  (B-l) 

5.  Calculation   o-f   New   Cycle   Engine   Expansion   Cylinder 
Preliminary  Design  Parameters: 

B,  L,  1 ,  a,  R,  NRpM,  max  Sp 

From  Appendix  A,   the  -following  preliminary  design  conditions 
Are    known: 

a.  Clearance  volume,  V_  =  0   :   L^,  =  1 ,_  +  a  -  s 

c  t»     r 

b.  Sleeve  clearance  =  0  :  B  =  bore  and  piston  diameter 

c.  Square  con-figuration   :   B  =  L 

d.  Sp  =  33.33  -ft /sec 

e.  VBDC  =  1.43333  -ft3 
*«  V900  =  0.80770  -ft3 


ill 


At  VBDC  Le  =  L  =  B: 


o 


VBDC  =  i-43333  "ft"  =  El*  Le  =  rrB-  <B-4) 

4         4 

B  =  L  =  C  4(1.45353)  31/3 

ir 

B  =  L  =  1.22204  -ft 

From  the  crank  radius  identity,  a  =  L/2: 

a  =  0.61102  ft 

From  equations  (B-l),  (B-4) ,  and  the  geometric  relationship, 
L^  =  l+a-sata  given  crank  angle  8: 

4V 

—   =  1   +  a  -  Cacose  +  ( 1  2-aji:sin2e)  X /23 
or  r 

ttB^ 

which  at  known  V,  B,  and  a  values  at  8  =  90°  simplifies  to: 

4 (.8077) 

-   =  lr+.61102-C.61102cos90°+(lr2-.611022sin290°) 1/2H 

tt(  1.22204)  2 

lr  -  .077615  =  (lr2-. 375546) 1/2 

lr2  -  .155251r  +  .006024  =  lr2  -  .575546 

1   =  2.4440  -ft 


ram  the  crank  ratio  identity,  Rc  =  lr/a: 
Rc  =  4 . 0 
From  equation  (B-5)  ,  NRp|y.  =  60  Sp/2L: 

60 ( 50 . 55 ) 


NRPM  = 


2(1.22204) 


NRpM  =  818.5 


From  equation  (B-2) ,  recognizing  that  maximum  Sp  occurs  at 
8  =  90° ,  270° : 

a2TTNRp|v)sin900  cos90° 


Sp     =   C  1+- 


max  60  (Rc2-sin290°) 1/: 
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Sp 


max 


=       a2ir_N 
60 


RPM 


(.61102)  (2rr)  (818.3) 
60 


=Rmax 


—      <5*3 


2.36    -ft/sec 
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FIGURE  B-1 :  Standard  reciprocating  cylinder  and  piston 
linkage  geometry 
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APPENDIX  C 

The  purpose  of  this  appendix  is  to  present  formulas  -for 

calculation   of   essential   thermal  properties   o-F   principle 

gases   as  a  function  of  mean  gas  temperature,   Tg ,   and   lean 

hydrocarbon-air   combustion   products,    C„H~    fuels,   as   a 

n  in        7 

function  of  Tg ,  pressure,  P,  and  equivalence  ratio,  $. 
1 .  Principle  gases 

TABLE  C-l 


Coefficients   for  calculation  of  C   and  >i  of  principle 
as  a  function  of  Tg<°R>  C63 

Valid  for  Tg  >  450° R 


gases 


Substance 


MW 

Y0 

104Y1 

108  Y2 

i°7JSa 

29 

,224 

.2617 

.  00 

2.  163 

28 

i— ,  ~^*^r 

.2478 

— 

.  00 1 2 

2.101 

28 

.229 

.  3306 

— 

.  2 1 60 

2.  053 

44 

.  141 

1.3128 

*-> 

.4867 

1 .  704 

18 

.479 

-.4483 

4 

.8951 

1.499 

air 

nitrogen  N^> 
carbon  monoxide  CO 
carbon  dioxide  C0o 
water  vapor  Ho0 


AMW  =  Atomic  Mass  Weight 


lbm/f t-sec 


=F  .645 

=  ^Qrg 


BTU/lbm-°R 


=  YQ  +  YxTg  +  Y2(Tg) 


(C-l) 


<C-2) 


2.  Lean  Hydrocarbon-Air  Combustion  products;  CnH^n  fuels  C73 
Valid  for  900° R  i  Tg  i  7200° R  and  $  &    1.0: 


1 bm/f t-sec 


1.470x10  7(Tg)"7 


l+.027$ 


(C-3) 
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FT  =  0.05  +  4.2Q-1)  -  6. 7(^-1  >*  (C-4) 

3T<Tg)  s  Figure  C-2 

Valid  for  900° R  £  Tg  £  3600° R,  P  =  1  atm,  and  §  =  . 75 
(see  Figure  C-2): 

C  =  .247  +  . 243xlO~4Tg  (C-5) 

PBTU/lbm-°R 
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FIGURE  C-1  :  Ratio  of  specific  heats,  0»  of  lean  hydrocarbon, 
air  combustion  products  as  a  function  of  ?,  Tg, 
and  $  ;  fuel  CnH2n  CO 
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FIGURE  C-2: 


Cp  of  lean  hydrocarbon-air  combustion  products 
as  a  function  of  P,  Tg,  and  $  J  fuel  CnH2n  C7J 
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